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
如有问题,请在评论区提问一起,一起交流进步。