均匀线阵方向图matlab程序代码

%线阵
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