%?MVDR?CAPON
clear;
i=sqrt（-1）;
j=i;
m=16;??%?m越大，分辨率越好
p=3;???%%?????????????
angle1=30;
angle2=-60;
angle3=45;
th=[angle1;angle2;angle3];
lma=2;
d=lma/2;
nn=1024;
SN1=10;
SN2=10;
SN3=10;
sn=[SN1;SN2;SN3];
degrad=pi/180;
%%construct?signals?and?noise
load?sig3
tt=1:nn
S=[s1（tt）.‘;s2（tt）.‘;s3（tt）.‘];
nr=randn（mnn）;
ni=randn（mnn）;
U=nr+j*ni;
Ps=S*S‘/nn;
ps=diag（Ps）;
refp=2*10.^（sn/10）;
tmp=sqrt（refp./ps）;
S2=diag（tmp）*S;
%%compute?covariance?matrix?and?eigenvalue?decomposition
tmp=-i*2*pi*d*sin（th‘*degrad）/lma;
tmp2=[0:m-1]‘;
a2=tmp2*tmp;
A=?exp（a2）;
X=A*S2+U;
Rxx=X*X‘/nn;
C=inv（Rxx）;
%%%spatial?spectrum
th2=[-90:1:90];
for?nnn=1:length（th2）
????tmp?=-i*2*pi*d*sin（th2（nnn）*degrad）/lma;
????tmp2?=?[0:m-1]‘;
????a2=tmp2*tmp;
????A2=exp（a2）;
????den=A2‘*C*A2;
????doa（nnn）=1/den;
end

%%%%%%?作出譜圖
semilogy（th2doa）;
axis（[-90?90?0.1?1e2]）;
xlabel（‘Angle（deg）‘）;
ylabel（‘Spectrum‘）;
grid;
hold?on

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????972??2010-08-19?10:20??陣列信號處理程序\CaponDOA.m

?????文件???????1922??2010-08-19?10:16??陣列信號處理程序\CaponDOA.asv

?????文件????????952??2010-08-19?10:20??陣列信號處理程序\music.m

?????文件???????1823??2010-08-18?20:53??陣列信號處理程序\singleSigSpec.m

?????文件???????1823??2010-08-18?20:52??陣列信號處理程序\singleSigSpec.asv

?????文件???????1190??2010-08-19?10:14??陣列信號處理程序\strongCorrelateMUSIC45and30degree.m

?????文件???????1191??2010-08-19?10:16??陣列信號處理程序\strongCorrelateMUSICminus60and30degree.m

?????文件???????1128??2010-08-19?10:11??陣列信號處理程序\strongCorrelateMUSICminus60and30degree.asv

?????目錄??????????0??2011-01-06?14:24??陣列信號處理程序

-----------?---------??----------?-----??----

????????????????11001????????????????????9