main.m

ESPRIT方法的两种实现

close all 
clear all  
format long 
N=128;
n=1:N;
k=20;
num=[0.200,0.213];
for i=1:k
    x=sqrt(20)*sin(2*pi*0.200*n)+sqrt(2)*sin(2*pi*0.213*n)+randn(1,N);
    Base1_LS(i,:)=Base_Esprit_1_Ls(x,12,50,num);
    Base1_TLS(i,:)=Base_Esprit_1_TLS(x,50,num);
    Base2_LS(i,:)=Base_Esprit_2_LS(x,12,50,num);
    Base2_TLS(i,:)=Base_Esprit_2_TLS(x,50,num);
end
mean(Base1_LS)
mean(Base1_TLS)
mean(Base2_LS)
mean(Base2_TLS)
var(Base1_LS)
var(Base1_TLS)
var(Base2_LS)
var(Base2_TLS)