%%ESPRIT角度和頻率聯合估計
%%BY胡德順??2015.1.10
%該方法比PM方法精度高很多

close?all
clear?all
clc

%%參數設置
k=3;
m=12;
p=3;
n=200000;
f=[35e549e54e6];
fc=5e6;
theta=[10*pi/18080*pi/18030*pi/180];
d=15;
t=1e-7;
fs=1e7;
c=3e8;
snr=10;
A=exp（-（j*2*pi*d/c）*（0:m-1）‘*（f.*sin（theta）））;
beta=exp（-j*2*pi*t*f）;
fai=diag（beta）;
S=sqrt（10^log（snr/10））*exp（j*2*pi*f‘*（1/fs）*（0:n-1））;
N1=randn（mn）;
N2=randn（mn）;
N3=randn（mn）;
X0=A*S+N1;
X1=A*fai*S+N2;
X2=A*fai*fai*S+N3;
X=[X0;X1;