阵列的方向图阵列输出的绝对值与来波方向图之间的关系称为天线的方向图。方向图一般有两类：一类是阵列输出的直接相加（不考虑信号及其来向），即静态方向图；另一类是带指向的方向图（考虑信号指向），信号指向是通过控制加权的相位来实现的。由信号模型可知，对于某一确定的元空间阵列，在忽略噪声的条件下，第个阵元的振幅为： (1)式中为来波的复振幅为第个阵元与参考点之间的延迟。设第个阵元的权值为，那么所有阵元加权的输出相加得到阵列输出为： (2)对上式取绝对值并归一化后可得到空间阵列的方向图： (3)如果式中，式（3）即是静态方向图。1 均匀线列阵假设均匀线阵间距为d，以左边的阵元为参考点，另假设入射方位角为，图1 均匀线列阵其中方位角表示与阵列法线方向的夹角，则有： (4)式(3)可以化简为： (5)式中，当时又可以进一步化简为： (6)可得均匀线阵静态方向图图： (7)当，的式(6)可以简化为： (8)于是可得指向为的阵列指向图： (9)MATLAB仿真图图2 指向0时，均匀线列阵的方向图图3 指向30时，均匀线列阵的方向图2. 均匀平面阵假设有一个的均匀面阵，其几何关系如图(3)，以阵列左上角的阵元为参考点，x轴上有n个间距为d的阵元，y轴上有m个间距为d的阵元。另假设信号的入射方位角为，俯仰角为，则信号入射到第个阵元上引起的与参考阵元间的时延为：图4 平面阵空间几何关系 (10)如果该阵元是竖直放置，则有： (11)当时的水平放置平面阵的方向图为： (12)竖直放置的平面阵方向图为： (12)即平面阵的方向图相当于合成行子阵与合成列子阵方向图的乘积。MATLAB仿真结果图5 平面阵水平和竖直放置时的指向性图3. 均匀圆阵若均匀线列阵由m个阵元组成，以均匀线列阵的圆心为参考点，圆的半径为r，时圆的指向性图为： (13)MATLAB仿真结果图6 均匀圆阵阵水平放置时的指向性图MATLAB仿真程序clear all;close all;clc;f=1000; %阵元发射信号频率c=1500; %声速lamta=c/f; %波长d=lamta/2; %阵元间距为半波长%均匀线阵%theta=-pi:pi/360:pi;phai=-pi:pi/360:pi; %方位角范围thetad=30*pi/180; %阵列指向betad=2*pi*d*sin(thetad)/lamta;m1=8; %阵元个数beta = 2*pi*d*sin(theta)/lamta;G11=sin(m1*(beta-betad)/2);G12=m1*sin(beta-betad)/2);G1=abs(G11./G12);GL=20*log(G1/max(G1);figureplot(theta*180/pi,GL);grid;axis(-180 180 -100 0);title(阵元M=8,thetad=30,均匀线阵方向图);xlabel(方位角/度);ylabel(G(theta)/dB);figurepolar(theta,G1);title(阵元M=8,thetad=30,均匀线阵方向图);%平面阵%x=16; y=16; %16*16面阵z=16;theta1=0:pi/360:pi; %方位角fai1=0:pi/360:pi; %俯仰角fai1,theta1=meshgrid(fai1,theta1);beta1=2*pi*d*cos(theta1).*cos(fai1)/lamta;beta2=2*pi*d*sin(theta1).*cos(fai1)/lamta;beta3=2*pi*d*sin(fai1)/lamta;Grow=0;for i1=1:1:x Grow1=exp(-i*(i1-1)*beta1); Grow=Grow+Grow1;endGcol=0;for k1=1:1:y Gcol1=exp(-i*(k1-1)*beta2); Gcol=Gcol+Gcol1;endGz=0;for j1=1:1:z Gz1=exp(-i*(j1-1)*beta3); Gz=Gz+Gz1;endGxy = Grow.*Gcol;Gxz = Grow.*Gz;figuremesh(theta1*180/pi,fai1*180/pi,10*log10(abs(Gxy/(x*y);grid on;axis(0 180 0 180 -60 0);xlabel(方位角theta/rad);ylabel(俯仰角phi/rad);zlabel(G(theta);title(水平放置平面阵方向图);figuremesh(theta1.*(180/pi),fai1.*(180/pi),10*log(abs(Gxz./(x*z);axis(0 180 90 180 -60 0);xlabel(方位角theta/rad);ylabel(俯仰角phi/rad);zlabel(G(theta);title(竖面放置平面阵方向图);%均匀圆阵%ff=100;lamta2=c/ff; %波长m3=16; %圆阵阵元数r=0.8090*lamta; %圆阵半径theta2=-pi:pi/360:pi; %方位角fai2=-pi:pi/360:pi; %俯仰角fai2,theta2=meshgrid(fai2,theta2);G=0;for mm=1:m3 G=G+exp(-i*(2*pi/lamta*r)*cos(2*pi*(mm-1)/m3-theta2).*cos(fai2);endfiguresurf(theta2*180/pi,fai2*180/pi,10*log10(abs(G./8);grid on;xlabel(方位角theta/rad);ylabel(俯仰角phi/rad);zlabel(G(theta);title(水平放置均匀圆阵方向图);