mvdr02.asv

窄带信号mvdr波束形成,这几天的学习结果

clc
clear all;
close all;
c=1500;
f0=6000;
fs=65536;  %采样频率
snap=1024;  %快拍数
lma=2;      %信号中心波长 
d=1;        %%天线阵元间距
M=8;        %------%基阵阵元数
p=4;        %信号个数
rad=pi/180;
theta=[0;50;110;80]*rad;
dd=cos(theta)*[0:M-1];
aa=[exp(2*pi*j*d*dd/lma)]';
n=1:snap;
v=[6000;6300;6500;6210];
ee=exp(j*v*2*pi*n);
snr=[0;10;5;15];
am=diag(sqrt(2)*10.^(snr/20));
noise=randn(M,snap)+j*randn(M,snap);
X=aa*am*ee;
R=X*X'/snap;
[E,S,V]=svd(R);           %计算特征值
Vs=E(:,1:p);          
Vn=E(:,p+1:M);              %小特征值对应的特征向量，即噪声子空间
r=inv(R);                   %求协方差矩阵的逆矩阵 

for i=0:180
    ap=[exp(2*pi*j*d*[0:M-1]*cos(i*rad)/lma)]';
    pm(i)=1/(ap'*Vn*Vn'*ap);
    pc
end