天线线列阵方向图

通信信号处理实验报告——阵列天线方向图的初步研究 11级通信（研） 刘晓娟一、实验原理：1、智能天线的基本概念：智能天线是一种阵列天线，它通过调节各阵元信号的加权幅度和相位来改变阵列的方向图形状，即自适应或以预制方式控制波束幅度、指向和零点位置，使波束总是指向期望方向，而零点指向干扰方向，实现波束随着用户走，从而提高天线的增益，节省发射功率。

智能天线系统主要由①天线阵列部分；②模/数或数/模转换部分；③波束形成网络部分组成。

本次实验着重讨论天线阵列部分。

2、智能天线的工作原理：智能天线的基本思想是：天线以多个高增益的动态窄波束分别跟踪多个期望信号，来自窄波束以外的信号被抑制。

3、方向图的概念：以入射角为横坐标，对应的智能天线输出增益为纵坐标所作的图称为方向图，智能天线的方向图有主瓣、副瓣等，相比其他天线的方向图，智能天线通常有较窄的主瓣，较灵活的主、副瓣大小、位置关系，和较大的天线增益。

与固定天线相比最大的区别是：不同的全职通常对应不同的方向图，我们可以通过改变权值来选择合适的方向图，即天线模式。

方向图一般分为两类：一类是静态方向图，即不考虑信号的方向，由阵列的输出直接相加得到；另一类是带指向的方向，这类方向图需要考虑信号的指向，通过控制加权相位来实现。

二、实验目的：1、设计一个均匀线阵，给出λ（波长），N （天线个数），d （阵元间距），画出方向图曲线，计算3dB 带宽。

2、通过控制变量法讨论λ，N ，d 对方向图曲线的影响。

3、分析旁瓣相对主瓣衰减的程度(即幅度比)。

三、实验内容：1、公式推导与整理：权矢量12(,,......)T N ωωωω=，本实验旨在讨论静态方向图，所以此处选择ω=(1,1,......1)T 。

信号源矢量(1)()[1,,...]j j N T a e e ββθ---=，2sin dπβθλ=，幅度方向图函数()()HF a θωθ==(1)1sin2sin 2Nj n n N eβββ--==∑=sin(sin /)sin(sin /)n d n d πθλπθλ。

阵列的方向图阵列输出的绝对值与来波方向图之间的关系称为天线的方向图。

方向图一般有两类：一类是阵列输出的直接相加（不考虑信号及其来向），即静态方向图；另一类是带指向的方向图（考虑信号指向），信号指向是通过控制加权的相位来实现的。

由信号模型可知，对于某一确定的m 元空间阵列，在忽略噪声的条件下，第l 个阵元的振幅为：l j l e g x ωτ-=0 (1)式中0g 为来波的复振幅l τ为第l 个阵元与参考点之间的延迟。

设第l 个阵元的权值为l w ，那么所有阵元加权的输出相加得到阵列输出为：∑=-=ml j l l e g w Y 100ωτ (2)对上式取绝对值并归一化后可得到空间阵列的方向图)(θG ：}max{)(00Y Y G =θ (3)如果式中1=l w ，式（3）即是静态方向图。

1． 均匀线列阵假设均匀线阵间距为d ，以左边的阵元为参考点，另假设入射方位角为θ，图1 均匀线列阵其中方位角表示与阵列法线方向的夹角，则有：d l cx c k l )1(1sin 1-==θτ (4)式(3)可以化简为：∑=--=ml l j l e g w Y 1)1(00β (5)式中λθπβ/sin 2=，当1=l w 时又可以进一步化简为：)2/sin()2/sin(2/)(00βββm m e mg Y l m j -= (6)可得均匀线阵静态方向图图：)2/sin()2/sin()(0ββθm m G =(7)当d l j l e w β)1(-=，λθπβdd d sin 2=的式(6)可以简化为：)2/)sin(()2/)(sin(2/)(00d d l m j m m e mg Y βββββ--=- (8)于是可得指向为d θ的阵列指向图：)2/)sin(()2/)(sin()(0d d m m G ββββθ--=(9)MATLAB 仿真图方位角/度G (θ)/d B2700阵元M=8,thetad=0,均匀线阵方向图图2 指向0°时，均匀线列阵的方向图方位角/度G (θ)/d B2700阵元M=8,thetad=30,均匀线阵方向图图3 指向30°时，均匀线列阵的方向图2. 均匀平面阵假设有一个n m ⨯的均匀面阵，其几何关系如图(3)，以阵列左上角的阵元为参考点，x 轴上有n 个间距为d 的阵元，y 轴上有m 个间距为d 的阵元。

阵列方向图与MATLAB 仿真1、线阵的方向图2()22cos(cos )R φψπφ=+-MATLAB 程序如下〔2元〕：clear;a=0:0.1:2*pi;y=sqrt(2+2*cos(pi-pi*cos(a)));polar(a,y); 图形如下：若阵元间距为半波长的M 个阵元的输出用方向向量权重11(,,)M j j M g eg e φφ⋅⋅⋅加以组合的话，阵列的方向图为 [(1)cos()]1()m Mj m m m R g e ψπφφ--==∑MATLAB 程序如下〔10个阵元〕：clear;f=3e10;lamda=(3e8)/f;beta=2.*pi/lamda;n=10;t=0:0.01:2*pi;d=lamda/4;W=beta.*d.*cos(t);z1=((n/2).*W)-n/2*beta* d;z2=((1/2).*W)-1/2*beta* d;F1=sin(z1)./(n.*sin(z2));iK1=abs(F1) ;polar(t,K1);方向图如下：2、圆阵方向图程序如下：clc;clear all；close all;M = 16; % 行阵元数k = 0.8090; % k = r/lambdaDOA_theta = 90; % 方位角DOA_fi = 0; % 俯仰角% 形成方位角为theta，俯仰角位fi的波束的权值m = [0 : M-1];w = exp(-j*2*pi*k*cos(2*pi*m'/M-DOA_theta*pi/180)*cos(DOA_fi*pi/180));% w = exp(-j*2*pi*k*(cos(2*pi*m'/M)*cos(DOA_theta*pi/180)*cos(DOA_fi*pi/180)+sin(2*pi*m'/M)*si n(DOA_fi*pi/180))); % 竖直放置% w = chebwin(M, 20) .* w; % 行加切比雪夫权% 绘制水平面放置的均匀圆阵的方向图theta = linspace(0,180,360);fi = linspace(0,90,180);for i_theta = 1 : length(theta)for i_fi = 1 : length(fi)a = exp(-j*2*pi*k*cos(2*pi*m'/M-theta(i_theta)*pi/180)*cos(fi(i_fi)*pi/180));%a=exp(-j*2*pi*k*(cos(2*pi*m'/M)*cos(theta(i_theta)*pi/180)*cos(fi(i_fi)*pi/180)+sin(2*pi*m'/ M)*sin(fi(i_fi)*pi/180))); % 竖直放置Y(i_theta,i_fi) = w'*a;endendY= abs(Y); Y = Y/max(max(Y));Y = 20*log10(Y);% Y = (Y+20) .* ((Y+20)>0) - 20; % 切图Z = Y + 20;Z = Z .* (Z > 0);Y = Z - 20;figure; mesh(fi, theta, Y); view([66, 33]);title('水平放置时的均匀圆阵方向图');% title('竖面放置时的均匀圆阵方向图'); % 竖直放置axis([0 90 0 180 -20 0]);xlabel('俯仰角/(\circ)'); ylabel('方位角/(\circ)'); zlabel('P/dB');figure; contour(fi, theta, Y);方向图如下：3、平面阵方向图：clc;clear all；close all;Row_N = 16; % 行阵元数Col_N = 16; % 列阵元数k = 0.5; % k = d/lambdaDOA_theta = 90; % 方位角DOA_fi = 0; % 俯仰角% 形成方位角为theta，俯仰角位fi的波束的权值Row_n = [0 : Row_N-1]; Col_n = [0 : Col_N-1];W_Row = exp(-j*2*pi*k*Row_n'*cos(DOA_theta*pi/180)*cos(DOA_fi*pi/180)); W_Col = exp(-j*2*pi*k*Col_n'*sin(DOA_theta*pi/180)*cos(DOA_fi*pi/180)); % W_Col = exp(-j*2*pi*k*Col_n'*sin(DOA_fi*pi/180)); % 竖直放置W_Row = chebwin(Row_N, 20) .* W_Row; % 行加切比雪夫权W_Col = chebwin(Col_N, 30) .* W_Col; % 列加切比雪夫权W = kron(W_Row, W_Col); % 合成的权值N*N x 1% 绘制水平面放置的平面阵的方向图theta = linspace(0,180,180);fi = linspace(0,90,90);for i_theta = 1 : length(theta)for i_fi = 1 : length(fi)row_temp = exp(-j*2*pi*k*Row_n'*cos(theta(i_theta)*pi/180)*cos(fi(i_fi)*pi/180)); % 行导向矢量N x 1col_temp = exp(-j*2*pi*k*Col_n'*sin(theta(i_theta)*pi/180)*cos(fi(i_fi)*pi/180)); % 列导向矢量N x 1% col_temp = exp(-j*2*pi*k*Col_n'*sin(fi(i_fi)*pi/180)); % 竖直放置Y(i_theta,i_fi) = W'*kron(row_temp, col_temp); % 合成的导向矢量N*N x 1 endendY= abs(Y); Y = Y/max(max(Y));Y = 20*log10(Y);Y = (Y+60) .* ((Y+60)>0) - 60; % 切图% Z = Y + 60;% Z = Z .* (Z > 0);% Y = Z - 60;figure; mesh(fi, theta, Y); view([66, 33]);title('水平面放置时的面阵方向图');axis([0 90 0 180 -60 0]);xlabel('俯仰角/(\circ)'); ylabel('方位角(\circ)'); zlabel('P/dB');figure; contour(fi, theta, Y);方向图如下：4、CAPON方法波束形成MATLAB程序如下〔阵元16，信号源3，快拍数1024〕：clear alli=sqrt(-1);j=i;M=16;%均匀线阵列数目P=3;%信号源数目f0=10;f1=50;f2=100;%信号频率nn=1024;%快拍数angle1=-15;angle2=15;angle3=30;%the signal angleth=[angle1;angle2;angle3]';SN1=10;SN2=10;SN3=10;%信噪比sn=[SN1;SN2;SN3];degrad=pi/180;tt=0:.001:1024;x0=exp(-j*2*pi*f0*tt);%3个信号x0、x1、x2x1=exp(-j*2*pi*f1*tt); %x2=exp(-j*2*pi*f2*tt); %t=1:nn;S=[x0(t);x1(t);x2(t)];nr=randn(M,nn);ni=randn(M,nn);u=nr+j*ni;%复高斯白噪声Ps=S*S'./nn;%信号能量ps=diag(Ps);refp=2*10.^(sn/10);tmp=sqrt(refp./ps);S2=diag(tmp)*S;%加入噪声tmp=-j*pi*sin(th*degrad);tmp2=[0:M-1]';a2=tmp2*tmp;A=exp(a2);X=A*S2+.1*u;%接收到的信号Rxx=X*X'./nn;%相关矩阵invRxx=inv(Rxx);%搜寻信号th2=[-90:90]';tmp=-j*pi*sin(th2'*degrad);tmp2=[0:M-1]';a2=tmp2*tmp;A2=exp(a2);den=A2'*invRxx*A2;doa=1./den;semilogy(th2,doa,'r');title('spectrum'); xlabel('angle'); ylabel('spectrum'); axis([-90 90 1e1 1e5]); grid;。

阵列天线方向图的MATLAB 实现课程名称：MATLAB程序设计与应用任课教师：周金柱班级：04091202姓名：黄文平学号：04091158成绩：阵列天线方向图的MATLAB 实现摘要：天线的方向性是指电磁场辐射在空间的分布规律，文章以阵列天线的方向性因子F（θ，φ）为主要研究对象来分析均匀和非均匀直线阵天线的方向性。

讨论了阵列天线方向图中主射方向和主瓣宽度随各参数变化的特点，借助M ATLAB绘制出天线方向性因子的二维和三维方向图，展示天线辐射场在空间的分布规律，表现辐射方向图的特点。

关键词：阵列天线;;方向图；MATLAB前言：天线是发射和接收电磁波的重要的无线电设备，没有天线也就没有无线电通信。

不同用途的天线要求其有不同的方向性，阵列天线以其较强的方向性和较高的增益在工程实际中被广泛应用。

因此，对阵列天线方向性分析在天线理论研究中占有重要地位。

阵列天线方向性主要由方向性因子F（θ，φ）表征，但F（θ，φ）在远区场是一组复杂的函数，如果对它的认识和分析仅停留在公式中各参数的讨论上，很难理解阵列天线辐射场的空间分布规律[ 1 ]。

MATLAB以其卓越的数值计算能力和强大的绘图功能，近年来被广泛应用在天线的分析和设计中。

借助MATLAB可以绘制出阵列天线的二维和三维方向图，直观地从方向图中看出主射方向和主瓣宽度随各参数的变化情况，加深对阵列天线辐射场分布规律的理解。

1 均匀直线阵方向图分析若天线阵中各个单元天线的类型和取向均相同，且以相等的间隔d 排列在一条直线上。

且各单元天线的电流振幅均为I，相位依次滞后同一数值琢，那么，这种天线阵称为均匀直线式天线阵，如图1 所示[ 2 ]:均匀直线阵归一化阵因子为[ 3 ]：Fn（θ，φ）是一个周期函数，所以除§= 0 时是阵因子的主瓣最大值外，§= ±2 mπ（m=1,2,...）都是主瓣最大值，这些重复的主瓣称为栅瓣，在实际应用中，通常希望出现一个主瓣，为避免出现栅瓣，必须把g限制在- 2π＜§＜2π范围内[ 4 ]，其中k=λ/2π，即波数，n 表示阵元数目。

目录一、基本概念 (1)1.1方向图基本概念 (1)1.2主瓣宽度 (2)1.2.1主瓣宽度基本概念及特性 (4)1.3旁瓣抑制 (4)一、基本概念1.1方向图基本概念天线的辐射电磁场在固定距离上随角坐标分布的图形，称为方向图。

用辐射场强表示的称为场强方向图，用功率密度表示的称之功率方向图，用相位表示的称为相位方向图。

天线方向图是空间立体图形，但是通常应用的是两个互相垂直的主平面內的方向图，称为平面方向图。

在线性天线中，由于地面影响较大，都采用垂直面和水平面作为主平面。

在面型天线中，则采用E平面和H平面作为两个主平面。

归一化方向图取最大值为一。

在方向图中，包含所需最大辐射方向的辐射波瓣叫天线主波瓣，也称天线波束。

主瓣之外的波瓣叫副瓣或旁瓣或边瓣，与主瓣相反方向上的旁瓣叫后瓣，见图1：全向天线水平波瓣和垂直波瓣图，其天线外形为圆柱型；图2：定向天线水平波瓣和垂直波瓣图，其天线外形为板状。

图1 全向天线波瓣示意图图2 定向天线波瓣示意图1.2主瓣宽度为了方便对各种天线的方向图特性进行比较，就需要规定一些特性参数。

主要包括：零功率波瓣宽度、半功率点波瓣宽度、旁瓣电平、前后比、方向系数等。

1.零功率波瓣宽度，指主瓣两侧场强值为0的两个方向之间的夹角，用2表示。

许多天线方向图的主瓣是关于最大辐射方向对称的，因此，只要确定零功率主瓣宽度的一半，再取其2倍即可求得零功率主瓣宽度,即2=2。

2. 半功率点波瓣宽度，指方向图主瓣两侧两个半功率点（即场强下降到最大值下降到0.707（或分贝值从最大值下降3dB处对应的两点）之间的夹角，又称为3dB波束宽度或主瓣宽度，记为。

对方向图对称天线，半功率主瓣宽度=2。

一般情况下，天线的E面和H面方向图主瓣宽度不同，分别记为、。

如不特殊说明，通常主瓣宽度是指半功率主瓣。

3. 副瓣电平，天线往往不止一个副瓣，而是有若干个。

仅靠主瓣的副瓣叫第一副瓣，依次为第二，第三、……副瓣，这些副瓣的峰值可能是不同的。

相控阵天线方向推导及仿真1、推导线阵天线方向图公式一个接收线阵，由等间距为d 的N 个各向同性单元组成，那么在θ方向，相邻单元接收信号的相位差为Ф=2πdλsinθ，线阵排列情况如图1所示。

图1 线阵排列示意图因为天线辐射方向图可以由天线上各种各样电流源辐射的单独贡献进行矢量叠加而得出，故各单元电压和为：E a =sin (ωt )+sin (ωt +ϕ)+sin (ωt +2ϕ)+⋯+sin⁡[ωt +(N −1)ϕ]将等式两边同时乘以2sin⁡(ϕ2)，根据积化和差、和差化积等相关数学公式，可得到如下公式：2sin (ϕ2)E a =cos (ωt −ϕ2)−cos (ωt +ϕ2)+cos (ωt +ϕ2)−cos (ωt −32ϕ)+⋯+cos (ωt +2N −32ϕ)−cos⁡(ωt +2N −12ϕ)整理得，2sin (ϕ2)E a =cos (ωt −ϕ2)−cos (ωt +2N−12ϕ)⁡⁡=2sin⁡(ωt +N −12ϕ)sin⁡(N2ϕ) 最终得到场强方向图，E a =sin⁡[ωt +(N −1)ϕ2⁄]sin⁡(Nϕ2⁄)sin⁡(ϕ2⁄)平方归一化后，得到辐射方向图（阵列因子）：|G a (θ)|=sin 2[Nπ(dλ)sinθ]N 2sin 2[π(dλ)sinθ]上式中，当(dλ)sinθ=0,±1,±2,···±n 时|G a (θ)|取得相等的最大值，但是我们只期望看到(dλ)sinθ=0的情况，取其他值产生的栅瓣是我们所不想见到的，为避免这种情况，特令d <λ。

前面的公式中认定主瓣指向为0°，当主瓣指向θ0方向时，则各向同性单元线阵的归一化辐射方向图为：G (θ)=sin 2[Nπ(dλ)(sinθ−sinθ0)]N 2sin 2[π(d λ)(sinθ−sinθ0)]此时，由于−2≤sin (θ)−sin (θ0)≤2，故防止产生栅瓣的条件为d <λ2⁄。

阵列天线方向图函数实验一、 实验目的1. 设计一个均匀线阵，给定d N d ,,,λθ画出方向图)(θF 函数图；2. 改变参数后，画出方向图)(θF 函数图，观察方向图)(θF 的变化并加以分析；3. 分析方向图)(θF 主瓣的衰减情况以及主瓣对第一旁瓣的衰减情况，确定dB3衰减对应的θ；二、 实验原理阵列输出的绝对值与来波方向之间的关系称为天线的方向图。

方向图一般有两类：一类是阵列输出的直接相加（不考虑信号及其来向），即静态方向图；另一类是带指向的方向图（考虑信号指向），当然信号的指向是通过控制加权的相位来实现的。

对于某一确定的M 元空间阵列，在忽略噪声的条件下，第k 个阵元的复振幅为),2,1(0M k e g x k j k Λ==-ωτ （2.1）式中：0g 为来波的复振幅，k τ为第k 个阵元与参考点之间的延迟。

设第k 个阵元的权值为k w ，那么所有阵元加权的输出得到的阵列的输出为),2,1(010M k e g w Y k j Mk k Λ==-=∑ωτ （2.2）对上式取绝对值并归一化后可得到空间阵列的方向图{}00max )(Y Y F =θ （2.3）如果),2,1(1M k w k Λ==式（2.3）即为静态方向图)(θF 。

下面考虑均匀线阵方向图。

假设均匀线阵的间距为d ，且以最左边的阵元为参考点（最左边的阵元位于原点），另假设信号入射方位角为θ，其中方位角表示与线阵法线方向的夹角，与参考点的波程差为θθτsin )1(1)sin (11d k cx ck -== （2.4）则阵列的输出为βθλπωτ)1(10sin )1(210100--=--=-=∑∑∑===k j Mk k d k jMk k j Mk k e g w eg w eg w Y k（2.5）式中：λθπβ/sin 2d =，λ为入射信号的波长。

当式（2.5）中),2,1(1M k w k Λ==时，式（2.5）可以进一步简化为)2/sin()2/sin(2)(00βββM M e Mg Y k M j == （2.6）可得均匀线阵的静态方向图，即)2/sin()2/sin()(0ββθM M F =（2.7）当式（2.5）中),2,1(,/sin 2,)1(M k d e w d d k j k d Λ===-λθπββ时，式（2.6）可简化为]2/)sin[(]2/)(sin[2)()1(00d d M j M M e Mg Y d ββββββ--=-= （2.7）于是可得到指向为d θ的阵列方向图，即]2/)sin[(]2/)(sin[)(d d M M F ββββθ--=（2.8）三、 实验过程1. 指向0=d θ静态方向图函数的实验1.1均匀线阵阵元个数N 对方向图函数)(θF 的影响sita=-pi/2:0.01:pi/2; lamda=0.03; d=lamda/2; n1=10; sita_d=0beta=2*pi*d*sin(sita)/lamda; beta_d=2*pi*d*sin(sita_d)/lamda; z11=(n1/2)*(beta-beta_d); z21=(1/2)*(beta-beta_d); f1=sin(z11)./(n1*sin(z21)); F1=abs(f1); figure(1);plot(sita,F1,'b'); hold on ; n2=20;beta=2*pi*d*sin(sita)/lamda;z12=(n2/2)*beta;z22=(1/2)*beta;f2=sin(z12)./(n2*sin(z22));F2=abs(f2);plot(sita,F2,'r');hold on;n3=30;beta=2*pi*d*sin(sita)/lamda;z13=(n3/2)*beta;z23=(1/2)*beta;f3=sin(z13)./(n3*sin(z23));F3=abs(f3);plot(sita,F3,'k')hold off;grid on;xlabel('theta/radian');ylabel('amplitude');title('¾ùÔÈÏßÕóÕóÔª¸öÊý¶Ô·½Ïòͼº¯Êý µÄÓ°Ïì'); legend('n1=10','n2=20','n3=30');分析：随着阵元数的增加，波束宽度变窄，分辨力提高。

阵列天线方向图及其MATLAB仿真一．实验目的1.了解阵列天线的波束形成原理写出方向图函数2.运用MATLAB仿真阵列天线的方向图曲线3.变换各参量观察曲线变化并分析参量间的关系二．实验原理1.阵列天线：阵列天线是一类由不少于两个天线单元规则或随机排列并通过适当激励获得预定辐射特性的特殊天线。

阵列天线的辐射电磁场是组成该天线阵各单元辐射场的总和—矢量和由于各单元的位置和馈电电流的振幅和相位均可以独立调整，这就使阵列天线具有各种不同的功能，这些功能是单个天线无法实现的。

^2.方向图原理：对于单元数很多的天线阵，用解析方法计算阵的总方向图相当繁杂。

假如一个多元天线阵能分解为几个相同的子阵，则可利用方向图相乘原理比较简单地求出天线阵的总方向图。

一个可分解的多元天线阵的方向图，等于子阵的方向图乘上以子阵为单元阵列天线天线阵的方向图。

这就是方向图相乘原理。

一个复杂的天线阵可考虑多次分解，即先分解成大的子阵，这些子阵再分解为较小的子阵，直至得到单元数很少的简单子阵为止，然后再利用方向图相乘原理求得阵的总方向图。

这种情况适应于单元是无方向性的条件，当单元以相同的取向排列并自身具有非均匀辐射的方向图时，则天线阵的总方向图应等于单元的方向图乘以阵的方向图。

三．源程序及相应的仿真图1.方向图随n变化的源程序clear;sita=-pi/2::pi/2;lamda=;]d=lamda/4;n1=20;beta=2*pi*d*sin(sita)/lamda;z11=(n1/2)*beta;z21=(1/2)*beta;f1=sin(z11)./(n1*sin(z21));F1=abs(f1);figure(1);plot(sita,F1,'b');hold on;n2=25;：beta=2*pi*d*sin(sita)/lamda;z12=(n2/2)*beta;z22=(1/2)*beta;f2=sin(z12)./(n2*sin(z22));F2=abs(f2);plot(sita,F2,'r');hold on;n3=30;beta=2*pi*d*sin(sita)/lamda;z13=(n3/2)*beta;z23=(1/2)*beta;>f3=sin(z13)./(n3*sin(z23));F3=abs(f3);plot(sita,F3,'k')hold off;grid on;xlabel('theta/radian');ylabel('amplitude');title('方向图与阵列个数的关系'); legend('n=20','n=25','n=30');·结果分析：随着阵列个数n的增加，方向图衰减越快，效果越好；2.方向图随lamda变化的源程序clear;sita=-pi/2::pi/2;n=20;d=;lamda1=;beta=2*pi*d*sin(sita)/lamda1;z11=(n/2)*beta;z21=(1/2)*beta;f1=sin(z11)./(n*sin(z21));~F1=abs(f1);%·½ÏòͼÇúÏßfigure(1);lamda2=;beta=2*pi*d*sin(sita)/lamda2;z12=(n/2)*beta;z22=(1/2)*beta;f2=sin(z12)./(n*sin(z22));F2=abs(f2);lamda3=;beta=2*pi*d*sin(sita)/lamda3;z13=(n/2)*beta;，z23=(1/2)*beta;f3=sin(z13)./(n*sin(z23));F3=abs(f3)plot(sita,F1,'b',sita,F2,'r',sita,F3,'k');grid on;xlabel('theta/radian');ylabel('amplitude');title('方向图与波长的关系');legend('lamda=','lamda=','lamda=');四．,随着波长lamda的增大，方向图衰减越慢，收敛性越五．结果分析：不是很好；3.方向图随d变化的源程序clear;sita=-pi/2::pi/2;n=20;lamda=;d1=;beta=2*pi*d1*sin(sita)/lamda;z11=(n/2)*beta;z21=(1/2)*beta;【f1=sin(z11)./(n*sin(z21));F1=abs(f1);%·½ÏòͼÇúÏßfigure(1);plot(sita,F1,'b');hold on;d2=;beta=2*pi*d2*sin(sita)/lamda;z12=(n/2)*beta;z22=(1/2)*beta;f2=sin(z12)./(n*sin(z22));F2=abs(f2);-plot(sita,F2,'r');hold on;d3=;beta=2*pi*d3*sin(sita)/lamda;z13=(n/2)*beta;z23=(1/2)*beta;f3=sin(z13)./(n*sin(z23));F3=abs(f3)plot(sita,F3,'k')hold off;grid on;xlabel('theta/radian');ylabel('amplitude');title('·½ÏòͼÓëÌìÏßÕóÁмä¸ôdµÄ¹Øϵ'); legend('d1=','d=','d=');结果分析；随着阵元之间间隔的增加，方向图衰减越快，主次瓣的差距越大，次瓣衰减越快，效果越好。

阵列天线方向图及其MATLAB 仿真1设计目的1.了解阵列天线的波束形成原理写出方向图函数2.运用MATLAB 仿真阵列天线的方向图曲线3.变换各参量观察曲线变化并分析参量间的关系2设计原理阵列天线：阵列天线是一类由不少于两个天线单元规则或随机排列并通过适当激励获得预定辐射特性的特殊天线。

阵列天线的辐射电磁场是组成该天线阵各单元辐射场的总和—矢量和由于各单元的位置和馈电电流的振幅和相位均可以独立调整，这就使阵列天线具有各种不同的功能，这些功能是单个天线无法实现的。

在本次设计中，讨论的是均匀直线阵天线。

均匀直线阵是等间距，各振源电流幅度相等，而相位依次递增或递减的直线阵。

均匀直线阵的方向图函数依据方向图乘积定理，等于元因子和阵因子的乘积。

二元阵辐射场：式中： 类似二元阵的分析，可以得到N 元均匀直线振的辐射场：令 ，可得到H 平面的归一化方向图函数，即阵因子的方向函数：式中：ζφθψ+=cos sin kd均匀直线阵最大值发生在0=ψ 处。

由此可以得出])[,(212121ζθθθϕθj jkr jkr m e r e r e F E E E E --+=+=12cos ),(21jkrm e F r E E -=ψϕθθζφθψ+=cos sin kd ∑-=+-=1)cos sin (),(N i kd ji jkrme erF E E ζϕθθϕθ2πθ=)2/sin()2/sin(1)(ψψψN N A =kdm ζϕ-=cos这里有两种情况最为重要。

1.边射阵，即最大辐射方向垂直于阵轴方向，此时 ，在垂直于阵轴的方向上，各元观察点没有波程差，所以各元电流不需要有相位差。

2.端射振，计最大辐射方向在阵轴方向上，此时0=mϕ或π，也就是说阵的各元电流沿阵轴方向依次超前或滞后kd 。

3设计过程本次设计的天线为14元均匀直线阵天线，天线的参数为：d=λ/2，N=14相位滞后的端射振天线。

基于MATLAB 可实现天线阵二维方向图和三维方向图的图形分析。