基于FDTD算法的波导色散特性仿真系统

资 源 简 介

本项目基于时域有限差分（FDTD）算法，在MATLAB环境下实现了对矩形波导及光纤等引导波结构的色散特性仿真与分析系统。其核心功能涵盖了计算电磁学的全流程，包括建立三维或二维Yee氏交错网格模型，并引入完美匹配层（PML）作为吸收边界条件以有效模拟无限大空间或边界损耗。通过在计算区域内施加高斯脉冲激励源，系统能够基于离散化的麦克斯韦旋转方程进行电场和磁场分量的时间步演进迭代，实时捕捉电磁波在波导内部的传播行为。该程序具备对观测点时域信号进行宽谱采集的能力，并利用快速傅里叶变换（FFT）将时域数据转换为频域

详 情 说 明

基于FDTD算法的波导色散特性MATLAB仿真研究

项目介绍

本项目是一个利用时域有限差分（FDTD）算法探究电磁波在引导波结构中传播特性的数值仿真系统。通过在MATLAB环境下构建二维计算模型，系统模拟了TE极化（包含电场分量Ey以及磁场分量Hx, Hz）在矩形介质波导中的演化过程。该工具专注于从时域数据中提取频域色散信息，能够精确计算传播常数、相速度及群速度，为微波器件与光波导的设计提供理论基础与仿真支撑。

功能特性

1. 全波数值演化：基于Maxwell旋转方程的离散化形式，实现电磁场在空域与时域的联合迭代更新。
2. 宽谱特性分析：通过高斯脉冲源激发，一次仿真即可获得宽频带范围内的响应数据。
3. 高效边界处理：在传播方向（Z轴）两端设置完美匹配层（PML），有效吸收传导波，消除人工边界反射。
4. 色散参数提取：利用相位差法处理双观测点信号，自动计算传播常数及归一化群速度曲线。
5. 动态交互可视化：仿真过程中实时显示电场强度分布图，并生成包含时域波形、透射频谱和色散关系的综合分析报表。

系统要求

1. 环境：MATLAB R2018a 或更高版本。
2. 工具箱：基础MATLAB功能即可运行，无需额外付费工具箱。
3. 硬件：建议4GB以上内存，以保证网格迭代时的实时显示流畅度。

核心功能与实现逻辑

系统的执行逻辑严格遵循计算电磁学的标准流程，具体步骤如下：

1. 参数初始化与网格划分 系统首先定义物理常数（如真空磁导率、介电常数及光速），并根据目标频率（THz量级）设定微米级的空间步长（dx, dz）。时间步长（dt）通过Courant稳定条件进行约束，确保数值计算的收敛性。

2. 几何结构与介质建模 在nx * nz的网格阵列中，通过修改介电常数矩阵定义波导结构。代码在中心区域建立了一个宽度受控的高折射率核心层（n=1.414），背景则默认为空气介质。

3. 吸收边界条件（PML）构建 为了模拟波向两侧无限延伸的物理环境，系统在Z轴两端应用了PML算法。通过设置四阶多项式分布的电导率（sigma_e）与磁导率（sigma_m），并依据波阻抗匹配原则进行关联，最大程度降低边界处的数值反射。

4. 激励源注入 系统采用硬源注入方式。激励源在时间上为高斯脉冲调制的正弦波，在空间上则被处理为与波导基本模态近似的正弦分布（Cosine Profile），以提高特定模态的激发效率。

5. 场分量迭代循环 这是仿真的核心模块，系统按时间步执行以下操作：

• 更新磁场Hx和Hz：利用当前步电场的空间差分计算磁场增量。
• 更新电场Ey：结合磁场的环流分布及介质的更新系数计算电场演化。
• 信号采集：实时记录两个预设观测点（obs_z1, obs_z2）处的Ey场强随时间的变化。

6. 后处理与色散提取 迭代完成后，系统对采集到的时域信号进行快速傅里叶变换（FFT）：

• 通过unwrap函数处理相邻观测点的相位差，获得不同频率对应的相位移动。
• 根据两点间的空间距离计算传播常数beta。
• 对beta随频率的变化率进行数值微分，求得群速度（v_g），能够分析结构色散对波速的影响。

算法实现细节

1. Yee氏交错网格：电场分量Ey位于网格边中心，磁场分量Hx, Hz则位于网格面中心，这种空间交错布局保证了旋度运算的中心差分精度。
2. 损耗系数同步：在更新系数Ce和Ch的计算中，充分考虑了介质电导率引起的损耗项，使得公式兼具普适性，不仅支持PML区域，也支持损耗介质的模拟。
3. 硬源处理技术：通过直接对特定网格位置的场值进行覆盖，确保了激励源的稳定输入。
4. 频域解析算法：采用相位平移法（Phase Shift Method）提取传播常数，避免了直接在空域求解特征方程的复杂性，能够直观反映宽带色散特性。