基于ESPRIT的DOA估计

clear all;close all;clc;



snapshot_number=500;             %快拍数
Ts=0.001;                        %采样周期
SNR=40;                          %信噪比


snr=15;%信噪比
tangle1=10;  %三个信号的发射角度
tangle2=35; %三个信号的入射角度
tangle3=60;

c1=kron（exp（-1i*pi*（0:9）*sin（tangle1*pi/180））‘exp（-1i*pi*（0:9）*sin（tangle1*pi/180））‘）;  %3个信号导向矢量     
c2=kron（exp（-1i*pi*（0:9）*sin（tangle2*pi/180））‘exp（-1i*pi*（0:9）*sin（tangle2*pi/180））‘）;       
c3=kron（exp（-1i*pi*（0:9）*sin（tangle2*pi/180））‘exp（-1i*pi*（0:9）*sin（tangle2*pi/180））‘）;

A=[c1c2 c3]; 
for r=1:50

s=5*（2*round（rand（31*snapshot_number））-1）;   %仿真信号
z1=A*s;
sigPower=mean（sum（abs（z1）.^2）/length（c1））;
noisePower=sigPower*10.^（-SNR/10）;
x=z1+randn（length（c1）snapshot_number）*sqrt（noisePower）;