阵列信号处理程序MUSIC MVDR 多信号分辨

用于雷达或移动通信或声纳中的阵列信号处理的MATLAB程序，包括多个信号的分辨力比较等

% 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