ESPRIT算法估计功率谱密度 matlab程序

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