clc;
clear all;
N=48;        %% N
k=1000;          %% sample number
thetas=10;%信號入射角度
thetai=[-30 30];%干擾入射角度
n=[0:N-1]';%n
sita=45*[-1:0.001:1];                                            %% scan vector
vs=exp(-j*pi*n*sin(thetas/180*pi));%信號陣列流型
vi=exp(-j*pi*n*sin(thetai/180*pi));%干擾流型
v=exp(-j*pi*n*sin(sita/180*pi));  %% 掃描流型 
f=15000;
t=[0:1:k-1]/1000;
snr=10;%信噪比
inr=10;%干噪比

xs=sqrt(10^(snr/10))*vs*exp(j*2*pi*f*t);%構(gòu)造有用信號
xi=sqrt(10^(inr/10)/2)*vi*[randn(length(thetai),k)+j*randn(length(thetai),k)];%構(gòu)造干擾信號(2個)
noise=[randn(N,k)+j*randn(N,k)]/sqrt(2);
% w=1/N*vs';
% figure(1),plot(sita,20*log10(abs(w*v)));
% axis([-90 90 -50 0]);
% hold on

%% mvdr 
NI=xi+noise;
Rni=NI*NI'/N;
wop1=Rni\vs/(vs'/Rni*vs);
B=abs(wop1'*v);
plot(sita,20*log10(B/max(B)),'r');
title('mvdr自適應(yīng)的波束圖(db)');xlabel('角度');ylabel('波束圖');
grid on
axis([-45 45 -50 0]);
hold off