N=2^16;
M=16;
f1=0.15;
f2=0.17;
f3=0.26;
SNR1=20;
SNR2=25;
SNR3=30;
A1=10^（SNR1/20）;
A2=10^（SNR2/20）;
A3=10^（SNR3/20）;
L=100;
f=[0.15?0.17?0.26];
est_f=zeros（13）;
est_error=zeros（13）;
xs=zeros（MN-M）;

for?l=1:L
????signal1=A1*exp（1i*（2*pi*f1*（0:N-1）+unifrnd（02*pi）））;
????signal2=A2*exp（1i*（2*pi*f2*（0:N-1）+unifrnd（02*pi）））;
????signal3=A3*exp（1i*（2*pi*f3*（0:N-1）+unifrnd（02*pi）））;
????noise=（randn（1N）+1i*randn（1N））/sqrt（2）;
????un=signal1+signal2+signal3+noise;
????
????for?k=1:N-M
????????xs（:k）=un（k+M-1:-1:k）.‘;
????end
????Rxx=xs（?:?1:end-1）*xs（?:?1:end-1）‘/（N-M-1）;
????Rxy=xs（?:?1:end-1）*xs（?:?2:end）‘/（N-M-1）;
????[UE]=svd（Rxx）;
????ev=diag（E）;
????emin=ev（end）;
????Z=[zeros（M-11）eye（M-1）;0zeros（1M-1）];
????Cxx=R