均匀阵列基于esprit算法的DOA估计

%基本ESPRIT算法第二种方法  最小二乘法
clear all;close all;clc;

c=3*10^8;
f=3*10^9;
%% 求得信号的波长
lamda=c/f;
%%　阵元的间距
d=lamda/2;
%%  （n-1）为子阵列的个数    即 阵元数
n=10;
%% 信号的数目
signal_number=3;
%% 三个信号的角度值
thita1=20;
thita2=30;
thita3=65;
%% 三个信号的中心频率
f1=40;
f2=20;
f3=70;
%% 在时域来说，是快拍数（一段时间内对阵列数据采样的个数）；在频域来说，是DFT的时间子段的个数。
snapshot=1:2000;
%% S是信号空间，由三个信号组成
S1=2.72*exp（1j*2*pi*f1*snapshot/length（snapshot））; 
S2=4.48*exp（1j*2*pi*f2*snapshot/length（snapshot））; 
S3=7.37*exp（1j*2*pi*f3*snapshot/length（snapshot））; 
S=[S1;S2;S3];
%% 子阵1
A1=exp（-1j*2*pi*d*[0:n-1]*sin（thita1*pi/180）/lamda）.‘;
A2=exp（-1j*2*pi*d*[0:n-1]*sin（thita2*pi/180）/lamda）.‘;
A3=exp（-1j*2*pi*d*[0:n-1]*sin（thita3*pi/180）/lamda）.‘;
A=[A1A2A3];

%%  噪声  假设为高斯白噪声，均值为零的
%N=zeros（102000）;
N= wgn（1020003）;
%% 求信噪比的  S1S2S3信噪比依次是10 20 30
s_power1=10*log（2.72^2/2）;
s_power2=10*log（4.48^2/2）;
s_power3