?
clc?
clear?all;?
close?all;??
N=16;?
lamda=1;?
d=lamda/2;?
theta0=0;?
theta=[-35*pi/18019*pi/18045*pi/180];?
INR=[403550];?
p=1;
%%%%%%%%%%%導向矢量%%%%%%%%%%%%%%%%%?
for?i=1:16?
atheta0（i1）=exp（j*（i-1）*2*pi/lamda*d*sin（theta0））;?
end?
for?i1=1:3?
p1（i1）=p*10^（INR（i1）/10）;?????
for?i=1:16?
atheta（ii1）=exp（j*（i-1）*2*pi/lamda*d*sin（theta（i1）））;?
end?
end?
R1=p*eye（N）;?
R=zeros（NN）;?
for?i1=1:3?
????tempt1=p1（i1）*atheta（:i1）*（atheta（:i1））‘;?
????R=R+tempt1;?
end?
R2=R;?
R=R+R1;?
%%%%%%%%%%%%%%%%%%%%%%%%%%?
E=eig（R）;%特征分解
E1=eig（R2）;?
[VD]=eig（R）;
k=1;?
for?theta1=-90:89?
?for?i=1:16?
atheta1（i1）=exp（j*（i-1）*2*pi/lamda*d*sin（theta1*pi/180））;?
end?
????Fq（k）=atheta0‘*atheta1（:1）;?
????k=k+1;?
end?
Fq1=abs（Fq）;?
l=1;?