【天线仿真】利用MATLAB实现多场景天线方向图可视化分析

1. 天线方向图到底是什么？从“手电筒”说起

如果你玩过手电筒，其实你已经理解了天线方向图的核心。把手电筒的光束想象成天线发射的无线电波，光束最亮、照得最远的方向，就是天线辐射最强的方向；而光束照不到的背后，就是天线辐射很弱甚至没有的方向。天线方向图，就是把这个三维立体的“光束形状”用图形画出来，告诉我们天线在空间的每个方向上“有多亮”。

在无线通信和雷达系统里，这个“光束形状”至关重要。比如，家里的Wi-Fi路由器，你希望它的信号（光束）能均匀覆盖整个客厅，而不是只对着墙壁猛照，这就需要一种类似“灯泡”的全向方向图。而像卫星电视的“小锅盖”，则需要把所有的能量像探照灯一样精准地聚焦到卫星的方向，这就需要尖锐的定向方向图。方向图直接决定了信号能不能传得远、覆盖得全、抗干扰能力强不强。

过去，工程师们靠公式计算和手工绘图，既抽象又繁琐。现在，有了MATLAB这样的工具，我们可以把复杂的电磁场计算变成直观的3D图形。你只需要写好描述天线特性的代码，比如天线单元的位置、电流大小，MATLAB就能帮你算出并画出那个“光束”，是胖是瘦，是圆是扁，一目了然。这就像给天线做了一次全面的“CT扫描”，设计好坏，一眼便知。

我刚开始接触时，也觉得那些辐射瓣、波束宽度、旁瓣电平的术语很唬人。但当你亲手用MATLAB跑出第一个方向图，看着那个立体图形旋转起来，你会瞬间明白：哦，原来这个参数调一调，主瓣就变窄了；那个参数改一改，旁边的杂散瓣（旁瓣）就压下去了。这种从代码到可视化的过程，是理解天线设计最直接、也最有成就感的方式。

2. 从零开始：你的第一个MATLAB方向图

别被理论吓到，我们从一段最简单的代码开始，目标是画出一个最基础的天线方向图。这里我们不直接用复杂的仿真引擎，而是用一个经典的数学模型——比如一个沿着Z轴放置的短偶极子天线——来生成数据并绘图。这会让你清晰地看到数据是如何流动，并最终变成图形的。

首先，我们需要在球坐标系中定义观察点的网格。天线辐射是三维的，所以我们用两个角度来定义空间中的任何一个方向：theta（俯仰角，从正Z轴往下量）和 phi（方位角，在XY平面内转圈）。在MATLAB里，我们可以这样创建一张覆盖整个空间的“角度网”：

% 定义角度网格
N_theta = 181; % theta采样点数，0到180度
N_phi = 361;   % phi采样点数，0到360度
theta = linspace(0, pi, N_theta); % 0到π，对应0到180度
phi = linspace(0, 2*pi, N_phi);   % 0到2π，对应0到360度
[Theta, Phi] = meshgrid(theta, phi); % 生成网格坐标

接下来，我们用一个简单的公式来计算这个短偶极子在每个方向上的辐射强度。对于沿Z轴的短偶极子，其功率方向图函数大约是 sin(theta)^2。注意，这个函数和 phi 无关，说明它在水平面是全向的。

% 计算辐射强度（假设为功率方向图）
radiation_pattern = (sin(Theta)).^2; % 这是场强幅度，功率则与它的平方成正比

现在，我们有了一个二维矩阵 radiation_pattern，它包含了在每一个 (theta, phi) 方向上的辐射强度值。怎么把它画成3D图呢？我们需要把球坐标 (r, theta, phi) 转换成MATLAB绘图能理解的直角坐标 (X, Y, Z)。转换公式是： X = R * sin(theta) * cos(phi) Y = R * sin(theta) * sin(phi) Z = R * cos(theta) 这里的 R 就是辐射强度的大小，它决定了图形在该方向上的“半径”。

% 将球坐标转换为直角坐标用于绘图
R = radiation_pattern; % 用辐射强度作为径向距离
X = R .* sin(Theta) .* cos(Phi);
Y = R .* sin(Theta) .* sin(Phi);
Z = R .* cos(Theta);

% 绘制3D方向图
figure('Position', [100, 100, 800, 600]); % 设置图形窗