%嘗試比較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）;???????????????????????????????????????