采用的是理论的归一化方向图公式
在代码里,theta0为天线波束指向方向,设置的步长为1,如果有30.5°等自行设置合适步长
%
clc;clear;close all;
c=3e8;f0=20e6;
lamda=c/f0;
d=lamda/2;
M=8;%阵元数
theta=-180:1:180;
theta0=0;
for i=1:M
a(i,:)=exp(1j*2*pi*(i-1)*d*(sind(theta)-sind(theta0))/lamda); %波束指向
end
F=sum(a);
[F_max,Index]=max(abs(F));
F_guiyi=abs(F)/F_max;
figure(1)
plot(theta,F_guiyi);title('8阵元线阵归一化方向图');xlabel('角度(°)');ylabel('幅值');
figure
polarplot(theta*pi/180,F_guiyi);title('8阵元线阵归一化方向图');
代码运行结果
在间距为半波长,8阵元,波束指向0时