mvdr02.m

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

%%1信号capon
clc
clear all;
close all;
c=1500;
f0=6000;
fs=65536;  %采样频率
snap=1024;  %快拍数
lma=c/f0;      %信号中心波长 
d=0.5*lma;        %%天线阵元间距
M=8;        %------%基阵阵元数
p=1;        %信号个数
rad=pi/180;
theta=[76]*rad;
dd=cos(theta)*[0:M-1];
aa=[exp(2*pi*j*d*dd/lma)]';
n=1:snap;
v=[6000];
ee=exp(j*v*2*pi*n);
snr=[-10];
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)=1/(ap'*Vn*Vn'*ap);
    pc(i+1)=inv(ap'*r*ap);
    pb(i+1)=ap'*R*ap;
end
pmm=10*log10(pm/max(pm));
pcc=10*log10(pc/max(pc));
pbb=10*log10(pb/max(pb))
figure(1);
plot(0:1:180,pmm);
figure(2);
plot(0:1:180,pcc);
figure(3);
plot(0:1:180,pbb);