天线方向图-阵元间距和波长的关系

clear
clc
N=8;%阵元数
d_lambda=[1/2,2];%d 和 lambda的关系
x=-0.5*pi:0.001*pi:0.5*pi;%角度范围
x0=0;%来波方向
p1=(1-exp(1j*2*pi*d_lambda(1)*N.*(sin(x)-sin(x0))))./(1-exp(1i*2*pi*d_lambda(1).*(sin(x)-sin(x0))));
p1=abs(p1);%P取模
p2=(1-exp(1j*2*pi*d_lambda(2)*N.*(sin(x)-sin(x0))))./(1-exp(1i*2*pi*d_lambda(2).*(sin(x)-sin(x0))));
p2=abs(p2);
for i=1:size(p1,2)
    if isnan(p1(i)) %替换为极限值
        p1(i)=N;
    end
end
for i=1:size(p2,2)
    if isnan(p2(i))
        p2(i)=N;
    end
end
p1_db=20*log10(p1./N);  %转化为dB
p2_db=20*log10(p2./N);
figure
subplot(1,2,1)
plot(x,p1_db);
xlim([-1.5 1.5])
xlabel('theta/radian')
ylabel('amplitude in dB')
title('d=λ/2')
grid on
subplot(1,2,2)
plot(x,p2_db);
xlabel('theta/radian')
ylabel('amplitude in dB')
title('d=2λ')
grid on

如有问题，请在评论区提问一起，一起交流进步。