圆形阵列的方向图matlab程序

%圆阵
clear all;
close all;
clc
N=16;%阵元个数
f=5000000000;%信号频率
c=300000000;%光速
lambda=c/f;%波长
d=0.5*lambda;%单元间距
theta=linspace（0360361）; %theta从0~360度扫描，共361个点
phi=linspace（09091）;%phi从0~90度扫描，共91个点
phi0=30;   %俯仰角
theta0=100; %预先设定的方向，方位角
Delta = 2*pi/N;%阵元之间的角度
R = d/Delta;%半径等于  弧长除以弧度
%R=8*lambda;
p = zeros（3N）;
k = 1;
for kk = 1:N
     p（:k） = [R*cos（kk*Delta）R*sin（kk*Delta）0]‘;
     k = k+1;
end   %阵元坐标
figure（1）;
plot3（p（1:）p（2:）p（3:）‘ko‘）;
hold on;
xlabel（‘\it x‘）;
ylabel（‘\it y‘）;
zlabel（‘\it z‘）;
k0=-2*pi/lambda*[sin（phi0*pi/180）*cos（theta0*pi/180）sin（phi0*pi/180）*sin（theta0*pi/180）cos（phi0*pi/180）].‘;
w0=exp（1i*k0.‘*p）;
for ii=1:length（theta）
    for jj=1:length（phi）
        k=-2*pi/lambda*[sin（phi（jj）*pi/180）*cos（theta（ii）*pi/180）sin（phi（