clear all
close all
N =10;
n = (-(N-1)/2:(N-1)/2)'
Vm = ones(N,1);
u=-1:1/1000:1;
INR=10^(10/10);
SNR=10^(20/10);
ui=0.15
us=0
Vi = exp(j*n*pi*ui)
Vs=exp(j*n*pi*us)

Sx=zeros(10,10)
%siganl_amp=sqrt(SNR);
%signal=signal_amp*exp(2*pi*j*randn)*(randn(10,1)+j*randn(10,1))/sqrt(2);
%x=signal*v
for i=1:100
X=sqrt(SNR)*exp(j*2*pi*randn)*(randn+j*randn)*Vi/sqrt(2)+exp(j*2*pi*randn)*(randn+j*randn)'*sqrt(INR)*Vs/sqrt(2)+(randn+j*randn)/sqrt(2)
Sx=(X*X')+Sx
end
Sx=Sx/100
%W=inv(Sx)*Vs/(Vs'*inv(Sx)*Vs)
%V=exp(j*n*pi*u)
u=zeros(101,1);
beam_db=zeros(101,1);
for i=1:201
    u(i)=(i-1)*0.01-1
V=exp(j*n*pi*u(i))
beam= 1/(V'*inv(Sx)*V);
beam_db(i)= 10*log10(abs(beam));

end
plot(u,beam_db)


%xlabel('{\itu}')
%ylabel('Beam2 pattern (dB)')
%title('\sigma_I^2 = 10dB, ui = 0.18')
%axis([-1 1 -30 10])
%line([0.15 0.15],[-30 5])
%grid