clear all
close all
A=1/10;
N =10;
n = (-(N-1)/2:(N-1)/2)';
Vm = ones(N,1);
u=1:100;
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)
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(SNR)*Vs/sqrt(2)+(randn+j*randn)/sqrt(2)
Sx=(X*X')+Sx
end
Sx=Sx/100
[U,S,V]=svd(Sx)
G=[];
for k=3:10
    G=[G,V(:,k)];
end
GG=G*G'
a=zeros(20,1)
for k=1:10
    for l=1:10
        a(10-k+l)=a(10-k+l)+GG(k,l);
    end
end
z=roots(a)
for k=1:19
  if abs(z(k))>1.02|abs(z(k))<0.98
    z(k)=0;
  end
end
z=sort(z);
anglez=angle(z)/2/pi;
absz=abs(z)
root1=abs(anglez(19)*2)
root2=abs(anglez(18)*2)