%圓陣
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（