天线辐射方向图二维三维可视化仿真系统

资 源 简 介

本项目是一个专为电磁场与微波技术初学者及工程技术人员设计的MATLAB仿真项目，旨在通过图形化手段直观展示简单天线的辐射特性。项目基于经典天线理论，重点实现了电赫兹偶极子（Hertzian Dipole）、半波偶极子（Half-wave Dipole）以及均匀直线阵列（Uniform Linear Array）等基础天线模型的方向图绘制。系统通过计算天线在空间不同角度下的辐射场强，利用MATLAB强大的绘图引擎，能够在二维平面上生成极坐标方向图和直角坐标方向图，用于分析E面和H面的辐射分布；同时，项目还具备三维可视化功能，能够构建出立体的空间辐射波瓣图，帮助用户从全方位理解天线的定向辐射能力。代码内部实现了自动归一化处理，能够清晰地展示主瓣、旁瓣和零点的位置。本项目不仅是一个实用的绘图工具，更是一个优秀的教学辅助资源，代码注释详尽，逻辑结构清晰，有助于用户深入理解天线增益、波束宽度、方向性系数等抽象概念，并掌握使用MATLAB进行科学计算和可视化的核心技能。

详 情 说 明

项目：简单天线方向图可视化教学演示系统

项目简述

本项目是一个基于MATLAB开发的仿真教学系统，旨在直观展示基础天线及阵列的辐射特性。系统专为电磁场与微波技术的初学者设计，通过交互式菜单引导用户选择不同的天线模型，利用MATLAB强大的数值计算与绘图能力，实时生成二维和三维的辐射方向图。该项目代码结构清晰，涵盖了从物理公式建模、数值网格生成到多维度数据可视化的完整流程，能够帮助用户深入理解波瓣宽度、旁瓣电平、主瓣方向等关键天线参数。

功能特性

系统主要具备以下核心功能：

1. 多模型支持：

* 电赫兹偶极子 (Hertzian Dipole)：模拟长度远小于波长的电流元辐射特性。 * 半波偶极子 (Half-wave Dipole)：模拟长度为半个波长的经典偶极子天线。 * 均匀直线阵列 (Uniform Linear Array)：支持用户自定义阵元数量、间距及扫描角，演示阵列天线的波束合成与扫描特性。

1. 多维度可视化：

* 2D 极坐标图：分别绘制E面（电场矢量平面）和H面（磁场矢量平面）的归一化方向图。 * 直角坐标功率谱：在直角坐标系下展示归一化功率（dB），便于观察波束宽度和旁瓣细节。 * 3D 立体波瓣图：构建三维空间辐射曲面，直观展示天线在全空间内的能量分布。

1. 自动参数分析：

* 自动计算并标注半功率波束宽度 (HPBW)。 * 自动识别最大辐射方向。 * 针对阵列天线，自动计算第一旁瓣电平 (SLL)。

系统要求

• 开发环境：MATLAB
• 工具箱依赖：系统主要依赖MATLAB基础绘图及数学函数，自动计算峰值部分使用了信号处理相关的基础逻辑（findpeaks函数），建议安装Signal Processing Toolbox以确保最佳兼容性，虽基础版本通常包含该函数的基本功能。
• 显示分辨率：建议在分辨率大于1280x720的屏幕上运行，以获得最佳的窗口布局体验。

详细代码实现与逻辑分析

本项目采用模块化设计，主程序通过无限循环配合图形化菜单实现持续交互，核心计算与绘图逻辑分离。以下是代码的具体实现细节分析：

1. 物理模型构建与计算逻辑

代码内部硬编码了工作频率为3GHz，并基于此计算波长和波数。针对三种不同的天线类型，分别实现了独立的计算函数：

• 电赫兹偶极子算法：

基于经典理论公式，场强与正弦函数 sin(theta) 成正比。代码生成网格后，直接计算场强绝对值，并进行归一化处理。为了防止对数计算出现负无穷，在计算dB值时对极小值进行了阈值限制（1e-6）。

• 半波偶极子算法：

实现了半波偶极子的标准方向函数公式，即分子为 cos(pi/2 * cos(theta))，分母为 sin(theta)。该模型展示了相比电赫兹偶极子更窄的主瓣宽度。

• 均匀直线阵列 (ULA) 算法：

* 参数配置：系统不仅是一个静态展示，还允许用户输入阵元数量(N)、以波长为单位的间距(d)以及波束扫描角。 * 相位控制：代码根据用户输入的扫描角，自动计算所需的单元间相位差 beta = -k * d * cos(theta0)，从而实现波束的电扫描。 * 阵因子计算：使用了通用的阵因子公式 sin(N*psi/2) / (N*sin(psi/2))。 * 奇异点处理：代码中包含关键的数值处理逻辑。由于阵因子分母在主极大方向会趋近于0，产生 0/0 的数值问题，代码通过洛必达法则逻辑，检测分母接近0的索引位置，强制将该点的阵因子值设为1，确保了数值计算的稳定性。 * 单元因子：为了突显阵列本身的特性，代码中将阵列单元假设为理想的各向同性源（Isotropic source），总场即为阵因子本身。

2. 多维度数据可视化逻辑

绘图核心函数 analyze_and_plot 负责将计算数据转化为可视图像，具体逻辑如下：

• 切面数据提取：

由于计算生成的是3D网格数据，绘图前先提取特定平面的切片数据。 * E面：提取方位角 phi = 0 对应的切片数据。 * H面：提取仰角 theta = 90度 对应的切片数据（垂直于偶极子轴线的平面）。

• 图形绘制布局：

使用 subplot(2, 2, x) 将窗口分为四个区域： 1. 左上 (2D 极坐标 E面)：使用 polarplot 绘制，并将0-pi的数据镜像拼接为0-2pi，形成完整的封闭曲线。 2. 右上 (直角坐标 dB图)：将角度转换为度数作为X轴，功率dB值作为Y轴。代码在此图中添加了辅助线，标记最大辐射方向和-3dB位置。 3. 左下 (2D 极坐标 H面)：展示水平面上的辐射分布，对于直立偶极子通常表现为圆形。 4. 右下 (3D 立体图)：核心逻辑是将球坐标系下的辐射强度转换为笛卡尔坐标系 (X, Y, Z)。转换公式中，坐标点离原点的距离由归一化场强 F_norm 决定，利用 surf 函数绘制出光滑的立体曲面，并配置了光照（Gouraud lighting）和视角，增强立体感。

3. 性能指标自动计算

系统具备简单的数据分析能力，主要针对E面数据进行分析：

• HPBW (半功率波束宽度)：通过搜索算法，从最大值位置向左、向右分别寻找功率下降至 -3dB 的角度位置，两者差值即为波束宽度。
• SLL (第一旁瓣电平)：代码使用了 findpeaks 函数寻找局部极大值。在排除主瓣（最大值）后，取第二大的峰值作为第一旁瓣电平。如果未检测到旁瓣（如简单偶极子），则该指标不显示。
• 结果展示：计算结果通过 annotation 文本框直接悬浮显示在绘图窗口的中心位置，方便用户读取。

使用方法

1. 启动MATLAB，将当前目录切换至项目所在文件夹。
2. 在命令窗口输入 main 并回车，或直接运行 main.m 文件。
3. 系统将弹出选择菜单：

* 选择 "电赫兹偶极子" 或 "半波偶极子" 将直接生成对应的可视化窗口。 * 选择 "均匀直线阵列" 将弹出一个对话框，需输入阵元数量（如8）、波长间距比（如0.5）和扫描角度（如90度为侧射，其他角度为扫描）。输入完成后点击确定生成图像。

1. 在生成的图窗中，可以旋转3D视图以从不同角度观察天线辐射波瓣。
2. 点击菜单中的 "退出系统" 即可结束程序。