MATLAB计算超表面的远场效果，多个图代替表征CST，HFSS仿真计算结果。 用仿真软件需要几个小时出结果，MATLAB可以几秒钟出结果，两者的结果是一样的。 可以计算三维远场，近场，theta，phi等等。 画图结果可以直接放到文献里面， 使用这种方法可以快速的出结果。 便于前期的仿真和后面的设计。

概述

本文介绍一套基于 MATLAB 实现的远场辐射特性仿真系统，专为编码超表面（Coding Metasurface）结构设计。该系统通过解析建模与数值计算相结合的方式，高效模拟超表面在特定工作频率下的反射远场方向图，可作为高频电磁仿真软件（如 CST、HFSS）结果的快速替代或预验证工具。其核心优势在于计算效率高、可视化能力强，并支持多维度（三维立体图、二维等高线图、方位/俯仰截面图、角度-幅度热力图）输出，便于工程分析与论文绘图。

系统功能架构

本系统主要包含以下功能模块：

1. 物理参数初始化
   用户设定工作频率（默认为 1 THz），系统自动计算波长、波数、单元周期等关键电磁参数，为后续建模提供基础。

1. 编码相位图定义
   通过预设的二维矩阵reflect_phi定义超表面每个单元的编码状态（0 或 1），进而映射为对应的反射相位（0 或 π 弧度）。该设计体现了“数字编码超表面”的基本思想：利用离散相位单元调控电磁波前。

1. 远场电场叠加计算
   采用惠更斯原理与远场近似，对每个单元的反射场在球坐标系中进行相干叠加。计算中考虑了：
   - 单元位置引起的传播相位延迟；
   - 单元中心作为辐射源点的几何模型；
   - 均匀幅度假设（所有单元反射幅度为 1）。

通过双重循环遍历所有单元，累加其在离散化的方位角（φ）与俯仰角（θ）网格上的复电场贡献，最终获得完整的远场复振幅分布E_total。

1. 多维可视化输出
   系统提供五类图形输出，全面表征远场特性：

• 图1：三维立体方向图
  将球坐标系下的电场模值转换为直角坐标系（Ex, Ey, Ez），绘制具有真实空间指向性的 3D 辐射瓣图，直观展示主瓣方向、旁瓣电平及波束宽度。

• 图2：二维归一化等高线图
  在 sinθcosφ–sinθsinφ 平面上绘制电场幅度的等高线填充图，便于观察波束在横向平面上的能量分布，适用于平面波束扫描分析。

• 图3：固定方位角的俯仰方向图
  用户可指定任意方位角 φ（如 180°），系统提取对应 φ 切片，绘制 θ 从 0° 到 90° 的一维反射幅度曲线，用于评估垂直面波束特性。

• 图7：固定俯仰角的方位方向图
  用户可指定任意俯仰角 θ（如 18°），系统提取对应 θ 切片，绘制 φ 从 0° 到 360° 的一维反射幅度曲线，用于评估水平面波束扫描能力。

• 图4 与 图5：角度-幅度热力图
  分别使用imagesc和pcolor绘制 θ–φ 平面上的电场幅度热力图，颜色深浅直观反映不同空间方向上的反射强度，适合快速识别主瓣位置与波束分裂现象。

1. 性能监控
   系统内置tic/toc计时机制，自动输出程序运行时间，便于评估计算效率与后续优化。

技术特点与适用场景

• 高效替代仿真：相比全波电磁仿真软件，本方法计算速度快（秒级），适合参数扫描、编码序列快速验证等场景。
• 物理模型清晰：基于远场叠加原理，模型透明，易于理解与扩展（如引入非均匀幅度、非周期结构等）。
• 绘图风格专业化：采用jet色图、无网格线、平滑着色（shading flat）等设置，使图形风格接近 CST/HFSS 输出，便于学术发表。
• 灵活可配置：用户仅需修改reflect_phi矩阵即可模拟任意 0/1 编码序列，支持大规模超表面（当前为 32×32 单元）。

总结

该 MATLAB 系统为编码超表面的远场性能评估提供了一套完整、高效、可视化的解决方案。它不仅可用于教学演示、论文绘图，还可作为智能超表面设计流程中的快速验证工具，显著提升研发迭代效率。未来可进一步扩展支持多频点、极化调控、非理想单元响应等高级功能，增强其工程实用性。

MATLAB计算超表面的远场效果，多个图代替表征CST，HFSS仿真计算结果。 用仿真软件需要几个小时出结果，MATLAB可以几秒钟出结果，两者的结果是一样的。 可以计算三维远场，近场，theta，phi等等。 画图结果可以直接放到文献里面， 使用这种方法可以快速的出结果。 便于前期的仿真和后面的设计。