MVDR波束形成MATLAB实现

clc
clear
%%设置初始参数
v=3e8;
f=2.4e9;
lamda=v/f;
d=0.5*lamda;
w=2*pi*f;
m=10;%阵元数
n=1000;%快拍数
snr=0;%0dB
inr=20;%20dB
thetas=90/180*pi;
thetai=40/180*pi;
%%生成导向向量a（theta）
theta=0:180;
theta=theta/180*pi;
tao=d*cos(theta)/v;
atheta = zeros(181,10);%阵列流形向量
for ii=1:181
    for jj=1:10
        atheta(ii,jj)=exp(1i*w*tao(ii)*(jj-1));
    end
end
atheta=atheta.';
%期望和入射信号导向向量产生
a_theta_s=zeros(