N=8;
n=[-(N-1)/2 : (N-1)/2]';
us=0;
u1=0.4,u2=-0.6;
Vs=exp(-j*pi*us*n);
V1=exp(-j*pi*u1*n);
V2=exp(-j*pi*u2*n);
p=20,q1=10,q2=20;
SNR=10^(p/10);
INR1=10^(q1/10);
INR2=10^(q2/10);
Rx=SNR*Vs*Vs'+INR1*V1*V1'+INR2*V2*V2'+eye(N);
w=inv(Rx)*Vs/(Vs'*inv(Rx)*Vs);
u=-1:1/1000:1;
beam=20*log10(w'*exp(-j*n*pi*u));
plot(u,beam);