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)
V=exp(j*n*pi*U)
Xk=[]
for i=1:1000
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)
Xk=cat(2,Xk,X);
end
P(:,:,1)=eye(10,10)/0.001;
W(:,1)=Vs/N;
u=0.5;
M(1)=1/(Vs'*P(:,:,1)*Vs);
for i=2:1000
    G(:,:,i)=(P(:,:,(i-1))*Xk(:,i))/(u+Xk(:,i)'*P(:,:,(i-1))*Xk(:,i));
    P(:,:,i)=(P(:,:,(i-1))-G(:,:,i)*Xk(:,i)'*P(:,:,(i-1)))/u;
    M(i)=1/(Vs'*P(:,:,i)*Vs);
    W(:,i)=(M(i)*(eye(10,10)-G(:,:,i)*Xk(:,i)')*W(:,(i-1)))/(u*M(i-1));
end
beam= W(:,1000)'*exp(j*n*pi*U);
beam= 20*log10(abs(beam));
plot(U,beam)
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