matlab仿真均匀线阵方向图(方位向和俯仰向)

matlab仿真均匀线阵方向图(方位向和俯仰向)

matlab仿真均匀线阵方向图(方位向和俯仰向)

本文使用Matlab对均匀线阵的方向图进行仿真，相关简单理论及公式见链接。由于链接里没有附上相关代码，我以自身理解分别对方位向和俯仰向进行了仿真，希望可以交流一下，发现问题。
思路主要是根据公式分别仿真天线增益和方位角、俯仰角的关系，然后二者相乘得到阵列方向图。

参考链接: https://blog.csdn.net/haoqibbwa_r/article/details/104498411

代码

j=sqrt(-1);
lamda=1;%波长
d=lamda/2;%阵元间隔
N=10;%阵元数
M=20;
dAz0=20*pi/180;%主瓣波束方位角
dEl0=30*pi/180;%主瓣波束俯仰角
dAz=(-90:89)*pi/180;
dEl=(-90:89)*pi/180;
F1=zeros(1,length(dAz));
for i=1:N
    F1=F1+exp(j*2*pi*(i-1)*d*(sin(dAz)-sin(dAz0))/lamda);
end
figure;plot(dAz*180/pi,10*log10(abs(F1)));xlabel('\theta(°)');ylabel('阵列增益/dB');title('均匀线阵方向图与方位角关系');
F2=zeros(1,length(dEl));
for i=1:M
    F2=F2+exp(j*2*pi*(i-1)*d*(sin(dEl)-sin(dEl0))/lamda);
end
figure;plot(dEl*180/pi,10*log10(abs(F2)));xlabel('\theta(°)');ylabel('阵列增益/dB');title('均匀线阵方向图与俯仰角关系');
F=zeros(length(dAz),length(dEl));
for i=1:length(dAz)
    for ii=1:length(dEl)
        F(i,ii)=F1(i)*F2(ii);
    end
end
figure;mesh(dEl*180/pi,dAz*180/pi,(abs(F)));ylabel('方位角(°)');xlabel('俯仰角(°)');title('均匀线阵方向图');

仿真结果如下图:

结果如图所示，我直接将阵列在方位向上和俯仰向上的天线增益相乘来表示阵列的方向图，不知是否可行，以及在平面阵列中，例如矩形阵列，天线增益在方位和俯仰上又是怎样表示？