clear all
close all
N =10;
n = (-(N-1)/2:(N-1)/2)';
Vm = ones(N,1);
u=1:100;
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);
Xs=zeros(10,1);
Xn=zeros(10,1)
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);
X=X/1.414;

Xm=exp(j*2*pi*randn)*(randn+j*randn)'*sqrt(INR)*Vs/sqrt(2);
Xs=cat(2,Xs,Xm)
Xm=sqrt(SNR)*exp(j*2*pi*randn)*(randn+j*randn)*Vi/sqrt(2)+(randn+j*randn)/sqrt(2);
Xn=cat(2,Xn,Xm)
Sx=(X*X')+Sx;
end
Sx=Sx/100;
WH=Vs'*inv(Sx)/(Vs'*inv(Sx)*Vs);
for j=1:100
SINR(j)=(WH*Vs*(Xs(:,(j+1))'*Xs(:,(j+1)))*Vs'*WH')/(WH*(Xn(:,(j+1))'*Xn(:,(j+1)))*WH');
end
m=0;
i=1;
for j=1:100
  m=SINR(j)+m;
  SINRNEW(i)=m/i;
i=i+1;
end
plot(u,SINRNEW)
%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