face_direction_pic.m

ULA1一个天线阵的方向图,可以用来估计信号的来波方向,face_direction_pic面阵的方向图 music,MN-music,capon,ESPRINT为四种进测量来波方向的算法 该程序包主

clear;close all;clc;
f=2e6;c=3e8;lamda=c/f;d=0.5*lamda;
N=64;            
%--------------------------产生坐标-----------------------------%
x=[0 d 2*d 3*d 4*d 5*d 6*d 7*d];y=x;
[X,Y]=meshgrid(x,y);
X=X';Y=Y';
X=reshape(X,1,N);Y=reshape(Y,1,N);
D=[X;Y];          
%--------------------------theta 方位,fai俯仰------------%
theta0=60*pi/180;fai0=60*pi/180;
%--------------------------产生导向矢量--------------------------%
k0=[cos(theta0)*cos(fai0)  sin(theta0)*cos(fai0)];
atheta0=exp(j*2*pi*k0*D/lamda).';
%-----------------------------------------------------------------%
p=0;
for theta=(-90:1:90)*pi/180;  
  p=p+1;
  q=0;
for fai=(0:1:90)*pi/180 ;
    q=q+1;
k=[cos(theta)*cos(fai)  sin(theta)*cos(fai)];
 atheta=exp(j*2*pi*k*D/lamda).' ;
 A(p,q)=atheta0'*atheta;
end 
end
C=max(max(abs(A)));
theta=-90:1:90;
fai=0:1:90;
mesh(fai,theta,abs(A)/C);
grid on;
figure
contour(fai,theta,abs(A)/C);
xlabel('俯仰')
ylabel('方位')
title('面阵方向图')
hold on;plot(fai0*180/pi,[-90:90],'r');
hold on;plot([0:90],theta0*180/pi,'r');