clc;
clear all;
c=3e8;
f=13e9;
N=10; %偶极子个数
lamda=c/f;
d=0.5*lamda; %0.25λ
xi=1*pi; %天线单元的相位差
k=2*pi/lamda;
phi=[0:2*pi/800:2*pi]; %方位角
framemax = 48;
for n=1:framemax
u=k*d*cos(phi)-xi/framemax*(n-1); %[0:pi]
F =abs(sin(N*(u/2))./sin(u/2)); %阵因子
end
pp=polar(phi,abs(F),'r')