clear?all;
clc;
close?all;

N=1000;
f1=0.1;
f2=0.25;
f3=0.27;
p=2*pi*rand（13）;
n=1:N;
s1（n）=sqrt（10^3）*exp（j*2*pi*f1*n+j*p（1））;
s2（n）=sqrt（10^3）*exp（j*2*pi*f2*n+j*p（2））;
s3（n）=sqrt（10^2.7）*exp（j*2*pi*f3*n+j*p（3））;
v=sqrt（1/2）*randn（1N）+j*sqrt（1/2）*randn（1N）;
x=s1+s2+s3+v;


for?i=1:N-3
????A（i:）=[x（i+3）x（i+2）x（i+1）x（i）];
end
B=A.‘*A‘.‘;
[V?D]=eig（B）;
D=eig（D）;