二维光子晶体PWE与FDTD能带及传输仿真平台

资 源 简 介

本项目旨在提供一个完整的二维光子晶体光学特性分析系统，结合了频域分析与时域动态模拟两种核心手段。 功能模块一采用平面波展开法（PWE），主要用于研究光子晶体在周期性势场下的本征状态。通过将空间变化的介电函数进行傅里叶级数展开，并结合布洛赫定理将麦克斯韦方程转化为矩阵特征值问题进行求解。该模块能够支持正方形、三角形等常见的晶格类型，精确计算不同偏振模式（TE模式和TM模式）下的色散关系曲线，即能带结构，进而识别出光子禁带的能量范围和位置。 功能模块二采用时域有限差分法（FDTD），重点解决光子晶体在实际应用

详 情 说 明

二维光子晶体能带结构及能量传输仿真平台

项目介绍

本项目是一款基于 MATLAB 开发的综合性光子力学分析工具，专注于二维光子晶体的物理特性研究。通过集成频域的平面波展开法（PWE）与时域的时域有限差分法（FDTD），本平台能够实现从材料能带设计到实际结构中电磁波传输行为的全流程模拟。

系统旨在帮助研究人员和学生直观地理解光子禁带的形成机理，并验证光子晶体波导、缺陷结构对光的调控作用。

---

功能特性

• 双模式仿真方案：融合了用于本征态分析的 PWE 算法和用于动态演化模拟的 FDTD 算法。
• 能带结构计算：支持计算光子晶体在特定布里渊区路径下的归一化频率分布，并自动识别光子禁带。
• 缺陷波导建模：通过在周期性结构中引入线缺陷（去掉一行介质柱），模拟光在受限空间内的传输。
• 完美匹配层（PML）：实现了基于电导率梯度的吸收边界，消除开放边界处的数值反射。
• 实时场图可视化：动态展示电磁波在结构中的传播过程，包含 Ez 场的瞬时分布。
• 频谱分析：通过快速傅里叶变换（FFT）计算透射能谱，量化频率响应特性。

---

运行环境与使用方法

1. 系统要求：安装有 MATLAB 2018b 或更高版本的个人电脑。
2. 基本操作：

* 直接运行主程序脚本。 * 程序将首先启动 PWE 计算循环，此时命令行会提示计算进度。 * 计算完成后，系统自动进入 FDTD 迭代阶段，并弹出一个包含三个子图的动态窗口。 * 仿真结束后，用户可以根据生成的关系曲线分析带隙位置与透射效率之间的对应关系。

---

系统实现逻辑与核心算法说明

整个仿真过程在单一主函数中完成，逻辑结构清晰，分为参数初始化、能带计算、时域仿真以及结果展示四个阶段。

#### 1. 物理参数与环境初始化 系统定义了晶格常数 $a=1mu m$，介质柱半径 $r=0.2a$。材料属性设定为高对比度背景（如硅与空气），相对介电常数差值构成了光子晶体的周期性势场。

#### 2. 平面波展开法（PWE）模块实现 该模块用于求解麦克斯韦方程组在倒阵空间中的特征值问题：

• 倒格矢构建：通过网格化生成平面波基向量组，决定了截断项数。
• 傅里叶变换：计算空间介电常数倒数 $eta(r)$ 的倒空间映射。对于圆形介质柱，程序利用一阶贝塞尔函数 $J_1$ 计算其形状因子（Shape Factor），精确描述介质分布。
• 路径扫描：沿着不可约布里渊区的特征路径（$Gamma to X to M to Gamma$）进行 $k$ 点扫描。
• 特征值求解：构造 TM 模式的矩阵算子，通过求解矩阵本征值获取归一化频率，进而绘制能带图。

#### 3. 时域有限差分法（FDTD）模块实现 该模块负责模拟波的实时传播：

• 网格构建：采用 $160 times 100$ 的空间网格，并将介质柱阵列映射到网格点的介电常数矩阵中。通过条件判断语句在 $y$ 轴中心位置留出无柱区域，形成线缺陷波导。
• Yee 算法迭代：交替更新 $Hx$, $Hy$ 磁场分量与 $Ez$ 电场分量。
• 边界条件：在仿真区域四周设置了 10 个网格宽度的 PML 吸收层，利用三次方规律变化的电导率分布来吸收向外传播的波。
• 激励源注入：采用高斯脉冲调制的正弦信号，其中心频率设定在光子带隙范围内，以测试波导的选择性传输特性。

#### 4. 可视化与数据后处理

• 动态展示：在 FDTD 循环中，每隔 20 个步长刷新一次电场分布图。
• 带隙渲染：程序会自动计算第一和第二能带之间的频率间隙，并在能带图中以灰色阴影区域标出。
• 透射谱计算：记录波导末端的监测点信号，结合源信号进行 FFT 变换。通过对比输出信号与输入信号的能量比值，生成归一化透射频谱图，直观反映波导对不同频率光的拦截或通过能力。

---

实现细节分析

• 数值稳定性：FDTD 步长 $dt$ 严格遵循库朗稳定条件（Courant Condition），确保波场更新不发生数值发散。
• 归一化处理：所有的频率均以 $omega a / 2pi c$ 的形式进行归一化，使得计算结果具有普适性，不依赖于绝对物理尺寸。
• 矩阵化运算：在 PWE 模块中大量使用 repmat 和矩阵乘法代替多重循环，显著提升了特征值问题的构造速度。
• 吸收性能：PML 的 sigma_max 经过优化，配合 10 层厚度，能将人工边界处的反射控制在较低水平，从而获得准确的稳态场分布。