为什么在matlab平面数组公式中实现出错

Question

我已尝试在 matlab 中实现如下所示的公式来创建如下所示的绘图。

但是当我使用下面的代码实现公式时，我得到了一个完全不同的图，如下所示。

我哪里做错了？

[theta,phi] = meshgrid(-pi:0.2:0.5*pi,-pi:0.2:pi);
f=10;
lambda=300/f;
k=(2*pi)/lambda;
N=10;
M=10;
dx=0.5*lambda;
dy=0.5*lambda;
ksi_x=k*dx.*sin(theta).*cos(phi);
ksi_y=k*dy.*sin(theta).*sin(phi);
AF=((1/M)*((sin(0.5*M*ksi_x))./(sin(0.5*ksi_x)))).*((1/N)*((sin(0.5*N*ksi_y))./(sin(0.5*ksi_y))));

surf(phi,theta,AF)
colorbar

Answer 1

您提供的图不是基于 (Θ, φ)，它们是 Ψ_x 和 Ψ_y。

只需将surf(phi,theta,AF) 替换为surf(ksi_y, ksi_x, AF) 即可。