1. 项目背景与核心问题
双缝干涉实验作为波动光学的经典范例，在光学教学和研究中具有不可替代的地位。传统解析方法在处理复杂边界条件时存在局限性，而时域有限差分法(FDTD)通过离散化麦克斯韦方程组，能够直观展现电磁场随时间演化的全过程。本项目聚焦FDTD方法在双缝干涉模拟中的两个关键技术：完美匹配层(PML)和理想磁导体(PMC)边界条件的实现。
PML作为吸收边界条件，其核心思想是通过引入各向异性损耗媒质，使入射电磁波在边界层内呈指数衰减。而PMC边界则模拟了磁场的完美反射特性，在实现对称性边界条件时尤为关键。Matlab作为工程计算的标准工具，其矩阵运算优势与FDTD的离散特性高度契合。
2. 仿真环境构建与参数设计
2.1 计算域离散化
采用Yee网格进行空间离散，电场和磁场分量在空间和时间上交错采样。对于典型可见光波段(λ=500nm)，网格尺寸Δx=Δy=λ/20=25nm可平衡计算精度与效率。时间步长根据Courant稳定性条件确定为Δt=Δx/(√2c)≈59.1as。
matlab复制% 网格参数示例
lambda = 500e-9; % 波长(m)
dx = lambda/20; % 空间步长
dt = dx/(sqrt(2)*3e8); % 时间步长
2.2 双缝结构建模
设计缝宽a=2μm，缝间距d=5μm的典型杨氏双缝。在计算域中心区域用PEC(理想电导体)定义不透明屏，留出两个矩形开口作为狭缝。为增强衍射效果，建议缝长至少为缝宽的5倍。
关键技巧：在网格对齐时，确保狭缝边缘与Yee网格的电场节点重合，避免阶梯近似带来的误差。
3. 边界条件实现细节
3.1 PML参数优化
采用Berenger提出的复坐标拉伸理论实现PML。对于10层PML，导电率分布采用多项式渐变：
σ(ρ) = σ_max*(ρ/d)^n
其中n=3~4可获得较好吸收效果，σ_max根据理论反射系数公式计算：
matlab复制R(θ) = exp[-2σ_max d cosθ / (n+1)ε0c]
实测表明，当PML厚度为λ/2时，对正入射波的反射率可低于-80dB。
3.2 PMC边界实现
在Matlab中通过磁场分量对称性实现PMC：
matlab复制Hx(1,:) = Hx(2,:); % 下边界
Hx(end,:) = Hx(end-1,:); % 上边界
这种处理相当于令边界处切向电场为零，法向磁场为零。
4. 激励源设置与场量更新
4.1 高斯脉冲源注入
采用软激励方式，在源平面注入宽带高斯脉冲：
matlab复制t0 = 20*dt; % 脉冲中心时间
tau = 10*dt; % 脉冲宽度
Ez_source = exp(-(t-t0)^2/tau^2);
脉冲频谱应覆盖感兴趣的光学频段，同时避免高频分量引起数值色散。
4.2 场量更新方程
依据麦克斯韦旋度方程离散化得到更新公式，以Ez分量为例：
matlab复制Ez(i,j) = Ceze(i,j)*Ez(i,j) + Cezh(i,j)*...
((Hy(i,j)-Hy(i-1,j))/dx - (Hx(i,j)-Hx(i,j-1))/dy);
其中系数矩阵Ceze、Cezh需根据介质属性预先计算。
5. 结果分析与可视化
5.1 干涉图样提取
在距离双缝D=50μm处设置观测线，记录稳态场分布。通过傅里叶变换将时域信号转换为频域：
matlab复制[Pxx,f] = periodogram(Ez_time_series,[],[],1/dt);
peak_idx = find(Pxx==max(Pxx)); % 找出主频
5.2 特征参数测量
- 条纹间距：Δy = λD/d
- 对比度：V = (Imax-Imin)/(Imax+Imin)
实测数据与理论值的偏差应小于5%，否则需检查边界反射或网格离散误差。
6. 常见问题排查指南
| 问题现象 | 可能原因 | 解决方案 |
|---|---|---|
| 仿真发散 | 时间步长过大 | 检查Courant条件，减小Δt |
| 边界反射明显 | PML参数不当 | 调整σ_max或增加PML层数 |
| 条纹不对称 | 网格未对齐 | 重新定位狭缝位置 |
| 高频振荡 | 激励源频谱过宽 | 增大高斯脉冲宽度tau |
7. 工程实践建议
- 内存优化：对于大型计算域，可采用分块计算策略，将场量数据分批写入硬盘。
- 并行加速：利用Matlab的parfor对空间网格循环进行并行化处理。
- 参数扫描：构建自动化脚本批量运行不同缝宽/间距组合，生成干涉图样数据库。
实际调试中发现，当PML与散射体距离过近时，会引入虚假反射。建议保持PML与最近散射体至少2λ间距。对于周期性结构，可采用Floquet边界条件替代PMC以提高计算效率。