DOA与频率联合估计算法的MATLAB仿真

%%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;