%線陣
clear?all;
close?all;
clc
N=16;%陣元個數
f=30000000;%信號頻率
c=30000000;
lambda=c/f;%波長計算公式
d=0.5*lambda;%單元間距
theta=linspace（0360361）;?
theta0=30;?%預先設定的方向
k?=?1;
p?=?zeros（3N）;
for?iii?=?0?:?N-1
????????p（:k）?=?[0iii*d0]‘;
????????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*[cos（theta0