波达方向DOA估计的Capon算法以及Music算法比较

%尝试比较music算法和capon算法的doa估计效果
clear all;close all;clc;
d=1;                                       %天线阵元间距
lamda=2;                               %信号中心波长
q1=pi/4;
q2=pi/3;
q3=pi/6;
q4=3*pi/4;                             %四个输入信号的波达方向
A1=[exp（-2*pi*j*d*[0:6]*cos（q1）/lamda）]‘;        %信号的方向向量    
A2=[exp（-2*pi*j*d*[0:6]*cos（q2）/lamda）]‘;
A3=[exp（-2*pi*j*d*[0:6]*cos（q3）/lamda）]‘;
A4=[exp（-2*pi*j*d*[0:6]*cos（q4）/lamda）]‘;
A=[A1A2A3A4];                                             %得出A矩阵（7行4列）
n=1:2000;                               %取2000个快拍
f1=0.015;
f2=0.05;
f3=0.02;
f4=0.035;                                %四个信号的频率
X=[1.3*cos（f1*n）;1*sin（f2*n）;1*sin（f3*n）;1*sin（f4*n）];      %输入信号矩阵
U=A*X;                                                                                %总的输入信号
C=cov（U‘）;                               %总输入信号的协方差矩阵                                       
[VD]=eig（C）;