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

c=3*10^8;
f=3*10^9;
%%?求得信號(hào)的波長(zhǎng)
lamda=c/f;
%%　陣元的間距
d=lamda/2;
%%??（n-1）為子陣列的個(gè)數(shù)????即?陣元數(shù)
n=10;
%%?信號(hào)的數(shù)目
signal_number=3;
%%?三個(gè)信號(hào)的角度值
thita1=20;
thita2=30;
thita3=65;
%%?三個(gè)信號(hào)的中心頻率
f1=40;
f2=20;
f3=70;
%%?在時(shí)域來(lái)說(shuō)，是快拍數(shù)（一段時(shí)間內(nèi)對(duì)陣列數(shù)據(jù)采樣的個(gè)數(shù)）；在頻域來(lái)說(shuō)，是DFT的時(shí)間子段的個(gè)數(shù)。
snapshot=1:2000;
%%?S是信號(hào)空間，由三個(gè)信號(hào)組成
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];

%%??噪聲??假設(shè)為高斯白噪聲，均值為零的
%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