AWG阵列波导光栅Matlab仿真：从原理到代码实现

简介阵列波导光栅AWG模型及其Matlab仿真项目是一份可直接用于课程设计与毕业设计的完整代码包面向通信工程、光电信息、电子信息等专业学生也适合希望入门光波导器件仿真的初学者。项目代码经过完整运行测试实现了光谱分析、模式求解、有效折射率计算等关键环节并提供交互式界面可在调整参数后即时观察输出光谱帮助理解器件设计与验证流程。压缩包共22个文件以17个m脚本为主另含1个mlx实时脚本、1个可安装应用及简要说明文档整体约1.07MB体积精简、结构清晰便于快速查阅和部署。目前已有173人参与学习配套说明文档提供使用指引代码支持二次修改与功能拓展对课程设计演示、毕业设计前期验证都有较高参考价值。1. 阵列波导光栅AWG的Matlab仿真到底要交付什么从一份课设任务说起毕设或课设做到“阵列波导光栅AWG模型及其在Matlab中的仿真”这道题最常见的翻车场景是解压拿到压缩包脚本能跑、图能出答辩时却说不出AWG为什么能把不同波长的光送到不同输出端口。这个题目真正要交付的不是一段能run的代码而是一个能解释清楚光学原理、能调参、能验证的Matlab仿真模型。AWG看着像集成光学的黑匣子实际上一旦拆成“衍射展开 相位累加 干涉聚焦”三个动作用高斯包络加相位叠加不到两百行代码就能复现八通道光谱响应。这篇文章适合光通信方向做毕设、集成光学课设以及想复现文献里AWG响应曲线但不想一上来就碰BPM的人。2. 把AWG光学原理拆成Matlab能算的表达式衍射级次、FSR与相位累加2.1 先想通AWG靠什么把波长“分开”AWG的结构可以拆成五段输入波导、输入平板波导也叫自由传播区、阵列波导、输出平板波导、输出波导。光从输入波导出来后在平板波导里自然衍射展开均匀耦合进一束长度递增的阵列波导。阵列波导每根的长度比相邻那根多一个固定的ΔL光走到阵列波导出口时波长不同、相位延迟就不同携带了“波长 × 长度差”的信息。最后经过输出平板波导的衍射叠加不同波长在输出面不同位置干涉相长形成空间上分开的谱线。一句话串起来AWG等效于一个衍射光栅加一个透镜阵列波导负责产生波长相关的相位差输出平板波导负责把这些相位差翻译成空间位置。理解这层逻辑后Matlab模型就可以绕过复杂的光场传播直接用“输入幅度包络 × 阵列波导相位 × 自由空间传播因子”叠加出来。很多同学一上来就扎进波动方程反而把这个最本质的关系丢了。2.2 仿真前必须锁定的四个参数λ₀、Δλ、FSR、m动笔写代码前先算出四个相互绑定的参数。中心波长λ₀通常取1550nmC波段光通信的基准。通道间隔Δλ常见0.8nm对应100GHz频率间隔。通道数Nch决定最后出几条谱线课设里4、8、16都有8通道最典型。关键在衍射级次m和自由谱范围FSR的关系。FSR是AWG能容纳的最大波长跨度超过这个范围光谱会周期性重复出现区分不出是哪个通道。简化关系是FSR λ₀² / (nc·ΔL)由于衍射级次m nc·ΔL / λ₀所以FSR ≈ λ₀ / m。同时阵列波导长度差ΔL由m决定ΔL m·λ₀ / nc。建模型前必须验算FSR是否能覆盖所有通道FSR要大于Nch×Δλ最好留出1.2倍以上余量否则边缘通道会被截断光谱出现折叠。这个校验是后面所有参数调优的前提很多课设代码“看着不对”就是在这埋下的雷。2.3 从光栅方程到阵列波导相位项k为什么必须带ncAWG的光栅方程写出来是ns·d·(sinθi sinθo) nc·ΔL m·λθi是输入衍射角θo是输出聚焦角。方程本身的含义是来自不同阵列波导的光在输出面某点干涉相长时光程差必须满足整数倍波长。落到Matlab里真正要算的只有两个量每个阵列波导入口的幅度以及它经过整段路径后的累计相位。相位项的写法是φ_j 2π·nc·(L0 j·ΔL) / λ。注意这里的λ是真空波长但波导内光走的等效光程必须乘nc所以波数k要写成k 2π·nc / λ。漏乘nc是最常见的错误后面会专门讲现象。先把参数计算写进脚本% 基本参数单位统一用米 lambda0 1.55e-6; % 中心波长 1550 nm nc 1.47; % 波导等效折射率 m 60; % 衍射级次先给一个合理值 dL m * lambda0 / nc; % 阵列波导长度差单位米 FSR lambda0 / m; % 自由谱范围单位米 Nch 8; % 通道数 dlambda 0.8e-9; % 通道间隔 100 GHz单位米 % 校验 FSR 是否覆盖所有通道 fprintf(dL %.4f um\n, dL * 1e6); fprintf(FSR %.2f nm\n, FSR * 1e9); fprintf(FSR / ch span %.2f\n, FSR / (Nch * dlambda));这段代码的作用是把物理参数先固化成一个自检结果。dL算出来约63.3μm如果是硅波导或氮化硅波导这个量级是合理的FSR约25.8nm覆盖8×0.86.4nm绰绰有余。FSR / ch span的输出如果小于1.2说明m选大了要往小调。参数说明m的初始值可以按F S R ≈ λ₀ / m估算想要FSR大于10nm就取m小于155。实际课设里m取40到80最常用太小会导致ΔL过短、相位差太小、光谱峰展不开太大则FSR不足以容纳通道数。3. 在Matlab里搭最小可运行模型两百行代码跑出8通道光谱3.1 模型选型为什么用高斯包络加菲涅尔叠加而不是BPM做AWG仿真常见做法有三条路光束传播法BPM、角谱法、标量衍射叠加模型。BPM是工程上的标准工具能算波导弯曲损耗、模式耦合等真实效应但课设里用它属于超配——要处理边界吸收层、模式初始化、步长收敛跑一次参数扫描动辄几分钟调参时特别磨人。角谱法需要补零、采样和多次傅里叶变换边界处理稍有不慎就出虚假反射。我一般建议课设和毕设直接用标量衍射叠加模型把输入波导在平板波导里的光场近似成高斯分布每根阵列波导只取它所在位置的幅度和累计相位最后在输出面做一次菲涅尔衍射叠加。这个模型保留了AWG所有宏观行为——通道间隔、FSR、旁瓣分布、串扰趋势足够解释原理、支持调参分析。要清醒知道它的边界不包含偏振效应、模场失配、波导弯曲损耗所以算出来的插损是相对值不是绝对值。但题目要的是“模型 仿真 解释”不是“与流片测试对标”这个模型完全够用。3.2 最小可运行代码参数、相位和输出面干涉核心代码分三段。第一段是几何与采样参数第二段是波长扫描主循环第三段从二维场里提取各通道光谱。全部贴出来能跑我习惯把这些写在一个脚本里方便调试。% 几何参数 M 81; % 阵列波导数量取奇数保证序号关于零对称 d 4e-6; % 阵列波导间距 4 um f 500e-6; % 平板波导焦距 500 um w0 30e-6; % 输入高斯束腰 30 um ns 1.0; % 平板区折射率简化模型 j (1:M) - (M 1) / 2; xj j * d; % 阵列波导入口横向坐标 A exp(-(xj / w0).^2); % 高斯包络幅度中心最强 % 输出面采样坐标范围要能覆盖 8 个通道加两侧旁瓣 xout -200e-6:0.5e-6:200e-6; % 波长扫描范围覆盖总带宽再留一两个 FSR 的余量 lambdaScan lambda0 (-4e-9:0.01e-9:4e-9); % 二维光谱场 spec2d zeros(length(xout), length(lambdaScan)); for ii 1:length(lambdaScan) lam lambdaScan(ii); k 2 * pi * nc / lam; % 阵列波导内波数注意带 nc phi k * (j * dL); % 每根阵列波导的累计相位 for ix 1:length(xout) r sqrt(f^2 (xout(ix) - xj).^2); % 阵列波导出口到输出面距离 E sum(A .* exp(1i * phi) .* exp(1i * 2 * pi * ns * r / lam) ./ sqrt(r)); spec2d(ix, ii) abs(E)^2; end end逻辑上分四层A是输入光耦合进每根阵列波导的幅度中心最强、两边衰减phi是波长相关的相位累加波长一变整组相位就变exp(i·2π·ns·r/λ)是光在平板区自由传播的相位积累1/sqrt(r)是二维柱面波的振幅衰减sum把所有阵列波导的贡献叠加起来就是干涉结果。参数说明M81是为了左右对称且波导数量足够填充高斯包络实际有效波导只有束腰宽度内的三四十根边缘波导幅度接近零对结果没影响但有计算开销。xout步长0.5μm要保证一个通道的谱峰位置至少有10个采样点后面提取通道光谱时才不会踩空。lambdaScan步长0.01nm可以看出谱线细节如果想快速验证改成0.05nm谱线形状仍在只是峰谷不够锐利。3.3 从二维场里提取8条通道光谱判断“物理上像不像”二维场spec2d是“输出面位置 × 波长”的强度矩阵要得到每条通道的光谱先算相邻通道在输出面上的空间间距。对AWG空间色散率近似D f·m / (ns·d)含义是单位波长的横向位移。通道位置xc (ch - (Nch1)/2) · D · Δλ正好把8个通道对称铺在零两侧。D f * m / (ns * d); % 空间色散率单位 m/m Pch zeros(Nch, length(lambdaScan)); % 每行一条通道的光谱 for ch 1:Nch xc (ch - (Nch 1) / 2) * D * dlambda; [~, idx] min(abs(xout - xc)); Pch(ch, :) spec2d(idx, :); end % 画归一化光谱单位 dB figure; plot(lambdaScan * 1e9, 10 * log10(Pch / max(Pch(:)) eps)); xlabel(波长 (nm)); ylabel(归一化功率 (dB)); grid on;这段代码的要点是不能用xout里恰好等于xc的点要用最近邻索引因为采样网格不一定正好落在通道中心。画出图后重点看三个特征峰值是否等间隔排列、相邻谱线的交叉点是否在同一功率水平、旁瓣是否低于主峰15dB以上。如果峰值间隔不是0.8nm说明D算错了或者dL不对如果旁瓣起伏不对称多半是输出面采样范围没覆盖完整干涉条纹。4. 参数怎么调才能让光谱立得住基准参数表与三步调参顺序4.1 可直接抄的基准参数表一套能出图的组合下面的参数组合是我反复调过的能稳定跑出八条间距均匀的谱线也方便你在此基础上做参数扫描分析。单位全部换算成米避免一半参数用微米一半用纳米。参数符号基准值说明中心波长λ₀1550 nmC波段基准通道间隔Δλ0.8 nm对应100GHz通道数Nch8课设最常见规模衍射级次m60决定FSR约25.8nm波导等效折射率nc1.47硅/氮化硅波导近似值平板区折射率ns1.0简化模型阵列波导数M81奇数保证对称阵列波导间距d4 μm直接决定通道间距平板焦距f500 μm决定输出面尺度高斯束腰w030 μm控制包络宽度这套参数下通道在输出面的间距约D·Δλ f·m·Δλ/(ns·d) ≈ 6μm8个通道总跨度42μm输出面采样范围±200μm足够显示完整的旁瓣结构。4.2 三步调参顺序先FSR再色散最后动M和w0调参最忌讳想到啥改啥我总结的顺序是固定三步每步只动一个变量跑完记录结果再动下一个。第一步锁FSR和衍射级次m。先验算FSR / (Nch×Δλ)是否大于1.2不够就把m调小。m从60调到50FSR从25.8nm变成31nm但ΔL也跟着变短相位差变小谱线会变宽所以m不是越大越好也不是越小越好。这个环节的目标只有一个保证边缘通道不塌陷。第二步定空间色散也就是让通道间距落在520μm。通道间距由D f·m/(ns·d)决定。间距太小时通道在输出面上挤在一起旁瓣互相干扰串扰变差间距太大则输出面范围要拉大同样采样点数下分辨率变差。最常见的手段是调整阵列波导间距d而不是改f因为f还会影响衍射叠加的相位曲率。d从4μm改成6μm通道间距就从6μm变成4μm谱峰更容易区分代价是输出面需要更大范围。第三步才是动M和w0。M决定旁瓣高度w0决定包络宽度。血泪经验是先把M固定在81、w0取2030μm跑出基准图确认通道间隔正确后再扫这两个参数否则两个变量一起动出了问题分不清是几何错了还是包络错了。4.3 用旁瓣电平判断模型够不够格-15dB与-25dB两道坎课设答辩最常被问的一句话是“你怎么知道模型是对的”。光凭谱线位置对不算完整旁瓣电平是一个直观的定量指标。判断标准我一般用两道坎旁瓣低于主峰15dB说明模型“能交差”主要物理行为都对低于25dB说明参数匹配得好可以写进报告做分析。影响旁瓣的主要是阵列波导数M和束腰w0。M越大干涉项越多旁瓣被压得越低但M超过一定数量后改善变缓计算时间线性上涨。w0的作用相反束腰越窄包络越集中等效参与干涉的波导数变少主峰变宽旁瓣变大束腰太宽则边缘波导分到的能量过低动态范围被压缩。实际操作时先固定w0把M从41加到161记录旁瓣下降曲线找到拐点再反过来调w0旁瓣收敛到-25dB附近就算到位。5. AWG仿真高频翻车点排查现象、原因与对应修复5.1 谱线全是NaN或者乱飞单位混用这是环境问题还是模型问题先看控制台有没有NaN警告。最常见的原因是半个脚本用微米、半个脚本用米比如lambda0写成1.55e-6但dlambda写成了0.8而不是0.8e-9导致通道间隔比中心波长还大干涉相位一团乱。解决方法是全程强制SI单位制所有长度用米所有波长用米只有出图时才乘1e6或1e9转成μm、nm。我习惯在脚本头部加一行注释“单位m”并在每个参数后面标注物理含义防止中途改参数时顺手填错。5.2 峰的位置对但通道间隔错k漏了折射率现象是谱线数量对、形态对但相邻峰间隔不是0.8nm整体谱线还朝一个方向偏移。原因九成是阵列波导相位里k用了2π/λ漏乘nc。阵列波导里光走的是介质中的相位积累光程必须乘等效折射率否则等效长度差ΔL对应的相位差整体偏小色散率也跟着变。解决办法是把两个波数分开写。阵列波导内用k_wg 2·pi·nc/lam平板区用k_slab 2·pi·ns/lam命名上区分开不要共用一个k。这样即使ns和nc取成一样代码逻辑也清晰答辩时解释起来也顺畅。5.3 串扰忽高忽低输出面采样点不够现象是相邻通道串扰不稳定有的通道到-30dB有的只有-10dB。这是典型的“通道峰值被采样网格漏掉”。通道位置xc是连续的采样网格是离散的最近邻索引可能偏离真实峰值半个网格。通道间距只有6μm时如果dx取2μm峰值定位误差能达到1μm串扰计算全部失真。把dx降到通道间距的十分之一以下也就是0.5μm左右问题基本消失。如果嫌计算慢可以分两遍第一遍用粗网格定位峰值附近第二遍只在峰值附近做加密采样两边结果对比确认一致后再用粗网格跑全谱。5.4 边缘通道功率塌陷、出现多组峰FSR容不下通道数现象是中间通道正常0通道和7通道功率明显偏低或者波长扫描范围里多出一组一模一样的峰。原因是FSR约等于甚至小于Nch×Δλ高阶衍射的光谱和低阶衍射的光谱在输出面重叠边缘通道的功率分散到两个位置。回到第2.2节的校验式FSR / (Nch×Δλ)如果接近1就要把m调小。m60时FSR25.8nm8通道只占6.4nm理论上富余但如果你改成16通道或者把Δλ改成1.6nm就必须重新验算。这也是我一直强调先把参数脚本跑起来再研究物理行为的原因边界条件错了后面所有分析都是空中楼阁。5.5 换一组参数就翻车没有做参数记录和回归验证现象很玄学昨天调的参数能出漂亮谱线今天把w0改了0.5μm整条谱线就变形了。其实不是玄学是没有记录参数、没有做回归验证。AWG的参数耦合关系很强M和w0同时微调就可能把旁瓣电平从-24dB推到-12dB。我建议每次跑完都在脚本末尾追加一行注释记下参数组和三个指标旁瓣、邻道串扰、边缘通道插损。格式类似下表可以直接抄进报告。Mw0 / μmd / μm旁瓣 / dB邻道串扰 / dB边缘插损 / dB81304-21.4-18.21.7121304-24.8-19.11.681204-18.9-16.52.3提示Matlab脚本中文注释乱码是老问题把.m文件用UTF-8编码保存再打开能避免大部分报错。导出图片用exportgraphics(gcf, awg_spectrum.eps, ContentType, vector)比print -depsc稳定字体也正常。6. 让AWG模型更有说服力二维光谱图、指标量化与文献对照6.1 用二维光谱图验证波长周期性与级次一维光谱图只看得到通道分配二维光谱图能直接看到AWG的级次周期性。把spec2d画成imagesc横轴是波长纵轴是输出面位置颜色是功率figure; imagesc(lambdaScan * 1e9, xout * 1e6, spec2d); set(gca, YDir, normal); xlabel(波长 (nm)); ylabel(输出面位置 (μm)); colorbar;图中会出现一组斜条纹每条斜条纹代表一个衍射级次。条纹的斜率就是空间色散率D相邻条纹的横轴间距就是FSR。答辩时把这个图放出来配合一句“周期性符合FSRλ₀/m的预期”比任何解释都直观。如果条纹不是斜线而是弯的说明f或ns的设置导致相位曲率异常回去检查输出面的距离r计算。6.2 把插损、串扰、3dB带宽量化成表光谱图画得再漂亮最终评审还是看数字。从Pch里可以一次算出三个常用指标% 各通道峰值功率与位置 [peakVal, peakIdx] max(Pch, [], 2); peakWl lambdaScan(peakIdx) * 1e9; % 相邻通道串扰在本通道峰值位置取邻道归一化功率 xtalk zeros(Nch, 1); for ch 1:Nch if ch 1 xtalk(ch) xtalk(ch) 10 * log10(Pch(ch, peakIdx(ch-1))); end if ch Nch xtalk(ch) xtalk(ch) 10 * log10(Pch(ch, peakIdx(ch1))); end end xtalk xtalk / 2; % 左右邻道取平均peakWl是每个通道峰值对应的波长拿它与理论值λ₀ (ch - (Nch1)/2)·Δλ对比误差小于0.01nm就说明通道分配正确。xtalk就是邻道串扰正常在-15dB以下。3dB带宽可以再加一行代码用插值算但不影响模型正确性作为加分项放附录即可。6.3 和文献响应做定性对照别纠结绝对精度最后一步是把你的光谱和经典AWG文献里的响应曲线做定性对比。文献里AWG通道响应的特征是三段式主峰、肩部、旁瓣你的模型复现出的形状应该类似只是旁瓣深度和肩部宽度有差异。差异来源很明确这个简化模型没有模场失配、波导侧壁粗糙度散射、弯曲损耗和偏振依赖所以插损绝对值和3dB带宽不能拿来做定量对标。对照的目的是验证趋势增大M旁瓣下降增大w0旁瓣上升、边缘通道插损改善增大d通道间距变大、FSR不变。把这些趋势做一两组扫描并画进报告你就把“模型自洽”这件事闭环了。我自己做这套模型时最后悔的是前期花了大量时间研究BPM的边界吸收层后来发现课设阶段一个衍射叠加模型足以把所有核心结论讲透。参数记录的习惯也是翻车两次后才养成的——现在每跑一组参数注释里自动带一条记录答辩前整理图表省了大量时间。希望帮到你。本文还有配套的精品资源点击获取