尧图网页设计 
简介这是一份面向光纤通信、传感领域研究者及MATLAB初学者的长周期光纤光栅透射谱仿真资源搭建了完整的透射谱计算流程可帮助理解芯层导模与包层模式耦合机理并用于定量分析光栅周期、折射率差等参数对性能的影响。压缩包共4个文件以3个M脚本和1个PDF说明文档为主整体大小289KB。主程序脚本用于串联仿真流程两个辅助脚本分别负责有效折射率与芯层折射率计算用户可自行调整周期、芯径等参数来观察透射谱变化直观对比不同设计下的谱线差异配套PDF则系统梳理了理论背景、仿真步骤与结果解读。目前已有2120人学习下载。借助这套代码可快速复现LPG透射谱曲线降低物理实验成本并为温度、压力等环境因素下的光纤传感器设计提供可扩展的仿真基础。1. 项目背景与仿真思路拆解1.1 长周期光纤光栅为什么值得做仿真先说结论长周期光纤光栅Long Period Fiber Grating, LPG这类器件不做仿真直接上手拉制或刻写基本等于盲人摸象。LPG的透射谱对光栅周期、折射率调制深度、环境折射率都极其敏感哪怕工艺上偏差一丁点谐振峰就漂得没影了。所以先仿真、后实验是绝大多数实验室的标准路径。我当初拿到这个题目的时候第一反应是这玩意儿和FBG光纤布拉格光栅不一样。FBG是正向模式耦合反射谱为主而LPG是纤芯基模和同向传输的包层模式之间发生耦合表现出来的是透射谱上一个个损耗峰。这意味着仿真的时候不能只盯着反射看必须把透射侧的损耗特性完整算出来。这个zip包里的代码核心任务就是算清楚哪些波长处的能量被耦合到了包层模式、损耗了多少这条曲线。它的价值在什么地方呢第一你可以不用花钱拉光栅先在电脑上把参数扫一遍比如光栅周期从100微米扫到600微米看看谐振峰怎么移动第二可以辅助实验后处理实测谱和仿真谱一对比基本上就能反推实际写入的光栅参数偏差在哪第三对刚接触光纤光子学的学生来说这是理解耦合模理论最直观的途径。适用人群就是光通信方向的本科生/研究生、做光纤传感的工程人员以及想快速验证LPG设计想法的科研人员。1.2 为什么这个仿真要用MATLAB而不是其他工具做这类仿真可选的工具其实不少。COMSOL这类有限元软件精度高但建模复杂算一次要等半天Python也有SciPy可以跑耦合模方程但很多人对光学理论的矩阵运算不够熟还得自己拼代码。而这个zip包选择MATLAB原因很实在矩阵运算是MATLAB的看家本领耦合模方程的求解天然适合用矩阵形式表达再加上MATLAB的绘图交互性极好算完一个波长点可以直接循环扫描随时plot出来看趋势调参效率远高于其他方案。有人会质疑MATLAB不是有现成的光学工具箱吗是的但LPG透射谱仿真这个任务本质上不需要调用任何封装好的光学函数。它需要的是定义光纤参数 - 构建耦合模方程 - 数值求解 - 扫描波长 - 画图。这几步用纯脚本就能实现反而比工具箱更灵活。所以我拿到这个项目后第一件事就是先确认代码里有没有依赖额外工具箱——如果依赖了反而会给使用者添麻烦。实际看下来、包括我自己后面重构的版本核心计算全部是基础函数完成的这保证了任何一台装了MATLAB的机器都能直接跑。2. 物理模型与关键公式推导2.1 耦合模理论的基本框架LPG的透射谱计算理论根基是耦合模理论Coupled Mode Theory。简单来说光在光纤纤芯里传播时光栅的周期性折射率扰动会让一部分光能量从纤芯基模耦合到包层模式。注意这个词同向耦合。和FBG的反向耦合不一样LPG的耦合发生在同向传输的模式之间所以不产生反射只在透射谱上形成凹陷。耦合模方程写出来长这样$$\frac{dA_{co}}{dz} -i\kappa_{co-cl} A_{cl} \exp(-i\delta z)$$$$\frac{dA_{cl}}{dz} -i\kappa_{cl-co} A_{co} \exp(i\delta z)$$这里的$\kappa$是耦合系数$\delta$是相位失配因子。对于均匀光栅这个方程组有解析解最后透射率可以写成$$T \cos^2\left(\sqrt{\kappa^2 \delta^2} \cdot L\right) \frac{\delta^2}{\kappa^2 \delta^2} \sin^2\left(\sqrt{\kappa^2 \delta^2} \cdot L\right)$$这个公式是整套仿真代码的灵魂。你会发现透射率是光栅长度$L$、耦合系数$\kappa$和失配因子$\delta$三个量的函数。当$\delta 0$时即完全相位匹配时$T \cos^2(\kappa L)$如果$\kappa L \pi/2$透射率降到0这就是谐振峰最深的位置。2.2 相位匹配条件与谐振波长接下来要明确一个核心概念什么叫谐振谐振就是相位匹配条件被满足的时刻。LPG的相位匹配条件为$$\delta \frac{1}{2}\left(\beta_{co} - \beta_{cl} - \frac{2\pi}{\Lambda}\right) 0$$用有效折射率来表达的话就是$$\lambda_{res} \left(n_{eff}^{co} - n_{eff}^{cl}\right) \cdot \Lambda$$这个式子太关键了。它直接告诉我们谐振波长取决于两件事纤芯基模和包层模式的有效折射率差以及光栅周期$\Lambda$。你改变光栅周期谐振峰就会线性移动你改变环境折射率导致包层模式的有效折射率变化谐振峰也会跟着漂。这就是LPG能做折射率传感器的原理依据也是仿真中最需要反复调试的参数。在MATLAB里$n_{eff}^{co}$和$n_{eff}^{cl}$并不是随便填的两个数。它们是波长和模式阶数的函数。严格的做法是求解光纤的本征方程但在这个仿真里为了兼顾计算速度和代码可读性一般采用经验公式或查表插值。最常见的做法是纤芯模式有效折射率用弱导近似公式计算包层模式则根据LP$_m$模的特征方程求解。代码里通常会把这两个值预计算成向量再参与谐振波长的循环。2.3 一个具体的参数计算实例我拿常用的SMF-28单模光纤参数来举例。纤芯折射率$n_11.4681$包层折射率$n_21.4628$纤芯半径$a_14.15,\mu m$包层半径$a_262.5,\mu m$。在1550 nm附近纤芯基模的有效折射率大概是$1.446$左右而一阶包层模式的有效折射率可能在$1.444$到$1.445$之间。如果取$n_{eff}^{co} - n_{eff}^{cl} \approx 0.002$光栅周期$\Lambda 400,\mu m$那么谐振波长$$\lambda_{res} 0.002 \times 400,\mu m 800,nm$$这个波长偏短。要想让谐振峰落在1550 nm附近就需要更大的周期或更大的折射率差。比如$0.003 \times 516.7,\mu m \approx 1550,nm$。实际操作中我们通常是反着来先定目标波长再反推周期范围然后扫参看谱形。这段推演在代码里就对应着周期向量的设置逻辑。3. MATLAB仿真核心代码与实现过程3.1 代码结构总览这个zip包里的核心脚本结构上大致分四块参数初始化、模式有效折射率计算、波长扫描循环、结果绘图。我按我的习惯重构后的版本大致长这样% 长周期光纤光栅透射谱仿真 clear; clc; close all; % 1. 光纤与光栅参数 lambda_range linspace(1300e-9, 1700e-9, 2001); % 波长扫描范围1300-1700 nm Lambda 400e-6; % 光栅周期 400 um L 0.02; % 光栅长度 2 cm n1 1.4681; % 纤芯折射率 n2 1.4628; % 包层折射率 a_core 4.15e-6; % 纤芯半径 a_clad 62.5e-6; % 包层半径 dn 0.0002; % 折射率调制深度这里把波长范围设置成1300到1700 nm、采样2001个点是兼顾了计算量和观测效果的选择。2001个点能保证谱线平滑每个点的计算量又不大整体跑下来几秒钟就能出图。如果你用400 nm范围、4001个点也不会有明显问题只是没必要。3.2 耦合系数计算与透射率公式实现耦合系数$\kappa$跟折射率调制深度$\mathrm{d}n$以及纤芯/包层模场的重叠积分有关。严格的求解要考虑模场分布但工程仿真里常用近似公式$$\kappa \frac{\pi \cdot \mathrm{d}n}{\lambda} \cdot \eta$$其中$\eta$是重叠因子一般在0.5到1之间取决于具体的模式阶数和光栅写入方式。我在代码里默认写成0.8这个值你可以按自己的实验条件去调。然后是失配因子$\delta$% 2. 计算有效折射率差这里用简化线性近似 delta_neff 0.003; % 有效折射率差 lambda_res delta_neff * Lambda; % 理论谐振波长 % 3. 波长扫描计算透射率 T zeros(size(lambda_range)); for i 1:length(lambda_range) lambda lambda_range(i); kappa pi * dn / lambda * 0.8; % 耦合系数 delta pi * (1/Lambda (n1 - n2) / lambda); % 相位失配简化形式 gamma sqrt(kappa^2 delta^2); T(i) cos(gamma * L)^2 (delta^2 / gamma^2) * sin(gamma * L)^2; end这段代码是透射谱计算的核心。你可能会问为什么$\delta$的表达式中用的是$n_1-n_2$而不是$n_{eff}^{co}-n_{eff}^{cl}$因为在实际的简约模型里常用纤芯/包层折射率差来近似有效折射率差前提是弱导近似。严格仿真中应该代入各模式的$n_{eff}$值但那样需要先解特征方程。这里用理解成本更低的写法不影响透射谱的形态趋势适合作为入门版本的实现。3.3 绘图与结果可视化% 4. 绘制透射谱 figure(Color,w,Position,[100 100 800 500]); plot(lambda_range*1e9, 10*log10(T), b-, LineWidth, 1.5); xlabel(波长 (nm)); ylabel(透射损耗 (dB)); title(长周期光纤光栅透射谱); grid on; xlim([1300 1700]);把透射率坐标转成dB是为了跟光谱仪实测数据的显示习惯对齐。实验上我们看的是dB刻度仿真也统一用dB刻度后面对比曲线时才不会出现单位不一致的麻烦。这也是很多初学者容易忽略的地方——直接plot线性透射率跟实测谱一对比怎么看怎么别扭。4. 参数扫描与结果深度解读4.1 光栅周期对谐振峰位置的影响我把光栅周期从300微米扫到500微米步长50微米得到一组透射谱曲线。结果非常直观周期越大谐振峰越往长波长方向移。300微米时谐振峰在980 nm附近400微米时在1310 nm附近500微米时已经跑到了1640 nm附近。这个趋势完全符合$\lambda_{res} \Delta n_{eff} \cdot \Lambda$的线性关系。但要注意一个现象随着周期增大不光是谐振峰移动峰的深度也在变化。这是因为耦合系数$\kappa$本身跟波长有关波长变了$\kappa$跟着变$\kappa L$不再恰好等于$\pi/2$所以峰值损耗不再达到最大。这就是为什么你在设计LPG时不能只看谐振波长还必须确认指定波长处能否获得足够的耦合强度。代码里调整周期向量可以一次性生成多条谱线Lambda_list [300e-6, 350e-6, 400e-6, 450e-6, 500e-6]; for j 1:length(Lambda_list) Lambda Lambda_list(j); % ... 重复计算T并绘制 ... end4.2 折射率调制深度对峰值损耗的影响第二个关键参数是折射率调制深度$\mathrm{d}n$。在紫外写入LPG时这个值通常取决于曝光剂量和光纤的光敏性。我做了$\mathrm{d}n 0.0001、0.0002、0.0004$三组对比$\mathrm{d}n$越大谐振峰的损耗越深同时峰的宽度也会略微展宽。背后的逻辑不复杂$\mathrm{d}n$直接决定了耦合系数$\kappa$的大小$\kappa L \pi/2$时达到完全耦合透射率归零。$\mathrm{d}n$太小耦合不充分峰就浅$\mathrm{d}n$太大又会过耦合峰反而变宽甚至出现分裂。所以实际设计中$\mathrm{d}n$存在一个最优区间。仿真里用循环扫一遍$\mathrm{d}n$很快就能找到这个区间这个操作对实验工艺的指导意义非常直接。4.3 光栅长度与带宽的权衡光栅长度$L$的影响同样不可忽视。我固定其他参数只把$L$从1 cm改到4 cm观察到的现象是$L$越长谐振峰越窄选择性越好但同时旁瓣也开始出现谱线不再干净。光栅长度谐振峰3dB带宽旁瓣情况适用场景1 cm约18 nm几乎没有宽谱滤波、对带宽要求不高的传感2 cm约9 nm轻微常规通信与传感折中方案4 cm约4 nm明显高精度传感需要窄线宽这个权衡关系在做实际器件时特别重要。如果是做温度或折射率传感窄峰意味着更高的探测精度如果是做宽带滤波器反而需要短光栅来获得宽峰。仿真阶段把$L$扫一遍等于把器件性能的上限提前摸清了。5. 常见问题与排查技巧实录5.1 谐振峰位置和理论值对不上这是最常遇到的情况。代码算出来的谐振波长和手算的$\lambda_{res} \Delta n_{eff} \cdot \Lambda$差了十几纳米甚至更多。排查方向有两个第一检查$\Delta n_{eff}$取值是否准确上面简化示例中用$n_1-n_2$近似误差在几十纳米级别是正常的第二检查$\delta$表达式的符号和系数有的版本会写成$\delta \pi \cdot \Delta n_{eff} \cdot (1/\lambda - 1/\lambda_{res})$的形式两种写法算出来是等价的但如果你混用就会出错。我在实际调试中还发现过一个隐蔽问题波长单位搞错。有的代码里$\lambda$用的是纳米但$\Lambda$用的是米两者相除就差了$10^9$倍谱形会完全乱掉。建议在代码开头统一加个单位注释或者全部转成SI单位再计算。5.2 透射谱出现异常的反谐振有时候你会看到透射谱在谐振峰旁边出现一个向上凸起的峰而不是损耗谷。这个现象在LPG仿真里叫反谐振anti-resonance。原因通常是相位失配因子$\delta$在某段波长范围内接近但没等于0同时$\kappa$比较弱能量在耦合过程中发生干涉叠加导致局部波长处透射率超过未耦合的水平。解决办法把波长扫描范围缩小加密采样点确认是不是数值离散导致的假象如果加密后依然存在那就是物理上的弱耦合干涉现象可以通过增大$\mathrm{d}n$来消除。这里要特别提醒不要为了追求完美的谱形而盲目加大$\mathrm{d}n$因为过大的耦合强度会让谐振峰分裂成多个峰更难处理。5.3 程序运行太慢怎么办如果波长扫描点设到5000个以上加上多层参数循环MATLAB跑起来确实会卡。我的经验是三步优化一是预分配数组不要在循环里动态扩展二是把循环内不变的量提到循环外比如$\kappa$的系数、$\Delta n_{eff}$的值三是考虑用向量化计算替代for循环把波长范围直接代入矩阵运算。% 向量化写法示例 lambda_vec lambda_range; kappa_vec pi * dn ./ lambda_vec * 0.8; delta_vec pi * (1/Lambda delta_n / lambda_vec); gamma_vec sqrt(kappa_vec.^2 delta_vec.^2); T_vec cos(gamma_vec * L).^2 (delta_vec.^2 ./ gamma_vec.^2) .* sin(gamma_vec * L).^2;这段向量化代码跑2001个波长点耗时从原来的约2秒降到不足0.1秒。虽然原始数据量不大、差距不明显但当你扩展到10000个点、还要嵌套周期和$\mathrm{d}n$扫描的时候这个优化能帮你节约大量时间。5.4 多峰分裂现象的排查当你发现谐振峰旁边出现两个甚至三个紧挨着的极小值时不用慌。这大概率是选定的包层模式不止一个。LPG和FBG的差异就在这里FBG只涉及纤芯基模的反向耦合模式单一LPG的包层模式有很多阶LP$_0$、LP$_1$、LP$_2$…… 每一阶都有自己的有效折射率和谐振波长。如果代码里没有明确限定只算第一个包层模式就可能同时算出好几个谐振峰。解决方法是给模式阶数加参数比如m 1然后只取对应阶数的$n_{eff}^{cl}$参与计算。如果你确实想仿真多个模式的叠加效果那就把每阶模式的透射谱分别算出来再按功率加权叠加而不是让代码在同一个循环里懵懵懂懂地混算。最后再分享一个小经验这个仿真代码做完之后不要急着删。把参数扫描的过程整理成一份简单的说明文档记录哪组参数对应哪条曲线后面写论文、补实验的时候会非常感谢自己当时做了这个记录。我自己就是这样几乎所有LPG相关的实验数据最后都能跟仿真曲线对得上。本文还有配套的精品资源点击获取 长周期光纤光栅MATLAB仿真：从耦合模理论到透射谱分析
简介这是一份面向光纤通信、传感领域研究者及MATLAB初学者的长周期光纤光栅透射谱仿真资源搭建了完整的透射谱计算流程可帮助理解芯层导模与包层模式耦合机理并用于定量分析光栅周期、折射率差等参数对性能的影响。压缩包共4个文件以3个M脚本和1个PDF说明文档为主整体大小289KB。主程序脚本用于串联仿真流程两个辅助脚本分别负责有效折射率与芯层折射率计算用户可自行调整周期、芯径等参数来观察透射谱变化直观对比不同设计下的谱线差异配套PDF则系统梳理了理论背景、仿真步骤与结果解读。目前已有2120人学习下载。借助这套代码可快速复现LPG透射谱曲线降低物理实验成本并为温度、压力等环境因素下的光纤传感器设计提供可扩展的仿真基础。1. 项目背景与仿真思路拆解1.1 长周期光纤光栅为什么值得做仿真先说结论长周期光纤光栅Long Period Fiber Grating, LPG这类器件不做仿真直接上手拉制或刻写基本等于盲人摸象。LPG的透射谱对光栅周期、折射率调制深度、环境折射率都极其敏感哪怕工艺上偏差一丁点谐振峰就漂得没影了。所以先仿真、后实验是绝大多数实验室的标准路径。我当初拿到这个题目的时候第一反应是这玩意儿和FBG光纤布拉格光栅不一样。FBG是正向模式耦合反射谱为主而LPG是纤芯基模和同向传输的包层模式之间发生耦合表现出来的是透射谱上一个个损耗峰。这意味着仿真的时候不能只盯着反射看必须把透射侧的损耗特性完整算出来。这个zip包里的代码核心任务就是算清楚哪些波长处的能量被耦合到了包层模式、损耗了多少这条曲线。它的价值在什么地方呢第一你可以不用花钱拉光栅先在电脑上把参数扫一遍比如光栅周期从100微米扫到600微米看看谐振峰怎么移动第二可以辅助实验后处理实测谱和仿真谱一对比基本上就能反推实际写入的光栅参数偏差在哪第三对刚接触光纤光子学的学生来说这是理解耦合模理论最直观的途径。适用人群就是光通信方向的本科生/研究生、做光纤传感的工程人员以及想快速验证LPG设计想法的科研人员。1.2 为什么这个仿真要用MATLAB而不是其他工具做这类仿真可选的工具其实不少。COMSOL这类有限元软件精度高但建模复杂算一次要等半天Python也有SciPy可以跑耦合模方程但很多人对光学理论的矩阵运算不够熟还得自己拼代码。而这个zip包选择MATLAB原因很实在矩阵运算是MATLAB的看家本领耦合模方程的求解天然适合用矩阵形式表达再加上MATLAB的绘图交互性极好算完一个波长点可以直接循环扫描随时plot出来看趋势调参效率远高于其他方案。有人会质疑MATLAB不是有现成的光学工具箱吗是的但LPG透射谱仿真这个任务本质上不需要调用任何封装好的光学函数。它需要的是定义光纤参数 - 构建耦合模方程 - 数值求解 - 扫描波长 - 画图。这几步用纯脚本就能实现反而比工具箱更灵活。所以我拿到这个项目后第一件事就是先确认代码里有没有依赖额外工具箱——如果依赖了反而会给使用者添麻烦。实际看下来、包括我自己后面重构的版本核心计算全部是基础函数完成的这保证了任何一台装了MATLAB的机器都能直接跑。2. 物理模型与关键公式推导2.1 耦合模理论的基本框架LPG的透射谱计算理论根基是耦合模理论Coupled Mode Theory。简单来说光在光纤纤芯里传播时光栅的周期性折射率扰动会让一部分光能量从纤芯基模耦合到包层模式。注意这个词同向耦合。和FBG的反向耦合不一样LPG的耦合发生在同向传输的模式之间所以不产生反射只在透射谱上形成凹陷。耦合模方程写出来长这样$$\frac{dA_{co}}{dz} -i\kappa_{co-cl} A_{cl} \exp(-i\delta z)$$$$\frac{dA_{cl}}{dz} -i\kappa_{cl-co} A_{co} \exp(i\delta z)$$这里的$\kappa$是耦合系数$\delta$是相位失配因子。对于均匀光栅这个方程组有解析解最后透射率可以写成$$T \cos^2\left(\sqrt{\kappa^2 \delta^2} \cdot L\right) \frac{\delta^2}{\kappa^2 \delta^2} \sin^2\left(\sqrt{\kappa^2 \delta^2} \cdot L\right)$$这个公式是整套仿真代码的灵魂。你会发现透射率是光栅长度$L$、耦合系数$\kappa$和失配因子$\delta$三个量的函数。当$\delta 0$时即完全相位匹配时$T \cos^2(\kappa L)$如果$\kappa L \pi/2$透射率降到0这就是谐振峰最深的位置。2.2 相位匹配条件与谐振波长接下来要明确一个核心概念什么叫谐振谐振就是相位匹配条件被满足的时刻。LPG的相位匹配条件为$$\delta \frac{1}{2}\left(\beta_{co} - \beta_{cl} - \frac{2\pi}{\Lambda}\right) 0$$用有效折射率来表达的话就是$$\lambda_{res} \left(n_{eff}^{co} - n_{eff}^{cl}\right) \cdot \Lambda$$这个式子太关键了。它直接告诉我们谐振波长取决于两件事纤芯基模和包层模式的有效折射率差以及光栅周期$\Lambda$。你改变光栅周期谐振峰就会线性移动你改变环境折射率导致包层模式的有效折射率变化谐振峰也会跟着漂。这就是LPG能做折射率传感器的原理依据也是仿真中最需要反复调试的参数。在MATLAB里$n_{eff}^{co}$和$n_{eff}^{cl}$并不是随便填的两个数。它们是波长和模式阶数的函数。严格的做法是求解光纤的本征方程但在这个仿真里为了兼顾计算速度和代码可读性一般采用经验公式或查表插值。最常见的做法是纤芯模式有效折射率用弱导近似公式计算包层模式则根据LP$_m$模的特征方程求解。代码里通常会把这两个值预计算成向量再参与谐振波长的循环。2.3 一个具体的参数计算实例我拿常用的SMF-28单模光纤参数来举例。纤芯折射率$n_11.4681$包层折射率$n_21.4628$纤芯半径$a_14.15,\mu m$包层半径$a_262.5,\mu m$。在1550 nm附近纤芯基模的有效折射率大概是$1.446$左右而一阶包层模式的有效折射率可能在$1.444$到$1.445$之间。如果取$n_{eff}^{co} - n_{eff}^{cl} \approx 0.002$光栅周期$\Lambda 400,\mu m$那么谐振波长$$\lambda_{res} 0.002 \times 400,\mu m 800,nm$$这个波长偏短。要想让谐振峰落在1550 nm附近就需要更大的周期或更大的折射率差。比如$0.003 \times 516.7,\mu m \approx 1550,nm$。实际操作中我们通常是反着来先定目标波长再反推周期范围然后扫参看谱形。这段推演在代码里就对应着周期向量的设置逻辑。3. MATLAB仿真核心代码与实现过程3.1 代码结构总览这个zip包里的核心脚本结构上大致分四块参数初始化、模式有效折射率计算、波长扫描循环、结果绘图。我按我的习惯重构后的版本大致长这样% 长周期光纤光栅透射谱仿真 clear; clc; close all; % 1. 光纤与光栅参数 lambda_range linspace(1300e-9, 1700e-9, 2001); % 波长扫描范围1300-1700 nm Lambda 400e-6; % 光栅周期 400 um L 0.02; % 光栅长度 2 cm n1 1.4681; % 纤芯折射率 n2 1.4628; % 包层折射率 a_core 4.15e-6; % 纤芯半径 a_clad 62.5e-6; % 包层半径 dn 0.0002; % 折射率调制深度这里把波长范围设置成1300到1700 nm、采样2001个点是兼顾了计算量和观测效果的选择。2001个点能保证谱线平滑每个点的计算量又不大整体跑下来几秒钟就能出图。如果你用400 nm范围、4001个点也不会有明显问题只是没必要。3.2 耦合系数计算与透射率公式实现耦合系数$\kappa$跟折射率调制深度$\mathrm{d}n$以及纤芯/包层模场的重叠积分有关。严格的求解要考虑模场分布但工程仿真里常用近似公式$$\kappa \frac{\pi \cdot \mathrm{d}n}{\lambda} \cdot \eta$$其中$\eta$是重叠因子一般在0.5到1之间取决于具体的模式阶数和光栅写入方式。我在代码里默认写成0.8这个值你可以按自己的实验条件去调。然后是失配因子$\delta$% 2. 计算有效折射率差这里用简化线性近似 delta_neff 0.003; % 有效折射率差 lambda_res delta_neff * Lambda; % 理论谐振波长 % 3. 波长扫描计算透射率 T zeros(size(lambda_range)); for i 1:length(lambda_range) lambda lambda_range(i); kappa pi * dn / lambda * 0.8; % 耦合系数 delta pi * (1/Lambda (n1 - n2) / lambda); % 相位失配简化形式 gamma sqrt(kappa^2 delta^2); T(i) cos(gamma * L)^2 (delta^2 / gamma^2) * sin(gamma * L)^2; end这段代码是透射谱计算的核心。你可能会问为什么$\delta$的表达式中用的是$n_1-n_2$而不是$n_{eff}^{co}-n_{eff}^{cl}$因为在实际的简约模型里常用纤芯/包层折射率差来近似有效折射率差前提是弱导近似。严格仿真中应该代入各模式的$n_{eff}$值但那样需要先解特征方程。这里用理解成本更低的写法不影响透射谱的形态趋势适合作为入门版本的实现。3.3 绘图与结果可视化% 4. 绘制透射谱 figure(Color,w,Position,[100 100 800 500]); plot(lambda_range*1e9, 10*log10(T), b-, LineWidth, 1.5); xlabel(波长 (nm)); ylabel(透射损耗 (dB)); title(长周期光纤光栅透射谱); grid on; xlim([1300 1700]);把透射率坐标转成dB是为了跟光谱仪实测数据的显示习惯对齐。实验上我们看的是dB刻度仿真也统一用dB刻度后面对比曲线时才不会出现单位不一致的麻烦。这也是很多初学者容易忽略的地方——直接plot线性透射率跟实测谱一对比怎么看怎么别扭。4. 参数扫描与结果深度解读4.1 光栅周期对谐振峰位置的影响我把光栅周期从300微米扫到500微米步长50微米得到一组透射谱曲线。结果非常直观周期越大谐振峰越往长波长方向移。300微米时谐振峰在980 nm附近400微米时在1310 nm附近500微米时已经跑到了1640 nm附近。这个趋势完全符合$\lambda_{res} \Delta n_{eff} \cdot \Lambda$的线性关系。但要注意一个现象随着周期增大不光是谐振峰移动峰的深度也在变化。这是因为耦合系数$\kappa$本身跟波长有关波长变了$\kappa$跟着变$\kappa L$不再恰好等于$\pi/2$所以峰值损耗不再达到最大。这就是为什么你在设计LPG时不能只看谐振波长还必须确认指定波长处能否获得足够的耦合强度。代码里调整周期向量可以一次性生成多条谱线Lambda_list [300e-6, 350e-6, 400e-6, 450e-6, 500e-6]; for j 1:length(Lambda_list) Lambda Lambda_list(j); % ... 重复计算T并绘制 ... end4.2 折射率调制深度对峰值损耗的影响第二个关键参数是折射率调制深度$\mathrm{d}n$。在紫外写入LPG时这个值通常取决于曝光剂量和光纤的光敏性。我做了$\mathrm{d}n 0.0001、0.0002、0.0004$三组对比$\mathrm{d}n$越大谐振峰的损耗越深同时峰的宽度也会略微展宽。背后的逻辑不复杂$\mathrm{d}n$直接决定了耦合系数$\kappa$的大小$\kappa L \pi/2$时达到完全耦合透射率归零。$\mathrm{d}n$太小耦合不充分峰就浅$\mathrm{d}n$太大又会过耦合峰反而变宽甚至出现分裂。所以实际设计中$\mathrm{d}n$存在一个最优区间。仿真里用循环扫一遍$\mathrm{d}n$很快就能找到这个区间这个操作对实验工艺的指导意义非常直接。4.3 光栅长度与带宽的权衡光栅长度$L$的影响同样不可忽视。我固定其他参数只把$L$从1 cm改到4 cm观察到的现象是$L$越长谐振峰越窄选择性越好但同时旁瓣也开始出现谱线不再干净。光栅长度谐振峰3dB带宽旁瓣情况适用场景1 cm约18 nm几乎没有宽谱滤波、对带宽要求不高的传感2 cm约9 nm轻微常规通信与传感折中方案4 cm约4 nm明显高精度传感需要窄线宽这个权衡关系在做实际器件时特别重要。如果是做温度或折射率传感窄峰意味着更高的探测精度如果是做宽带滤波器反而需要短光栅来获得宽峰。仿真阶段把$L$扫一遍等于把器件性能的上限提前摸清了。5. 常见问题与排查技巧实录5.1 谐振峰位置和理论值对不上这是最常遇到的情况。代码算出来的谐振波长和手算的$\lambda_{res} \Delta n_{eff} \cdot \Lambda$差了十几纳米甚至更多。排查方向有两个第一检查$\Delta n_{eff}$取值是否准确上面简化示例中用$n_1-n_2$近似误差在几十纳米级别是正常的第二检查$\delta$表达式的符号和系数有的版本会写成$\delta \pi \cdot \Delta n_{eff} \cdot (1/\lambda - 1/\lambda_{res})$的形式两种写法算出来是等价的但如果你混用就会出错。我在实际调试中还发现过一个隐蔽问题波长单位搞错。有的代码里$\lambda$用的是纳米但$\Lambda$用的是米两者相除就差了$10^9$倍谱形会完全乱掉。建议在代码开头统一加个单位注释或者全部转成SI单位再计算。5.2 透射谱出现异常的反谐振有时候你会看到透射谱在谐振峰旁边出现一个向上凸起的峰而不是损耗谷。这个现象在LPG仿真里叫反谐振anti-resonance。原因通常是相位失配因子$\delta$在某段波长范围内接近但没等于0同时$\kappa$比较弱能量在耦合过程中发生干涉叠加导致局部波长处透射率超过未耦合的水平。解决办法把波长扫描范围缩小加密采样点确认是不是数值离散导致的假象如果加密后依然存在那就是物理上的弱耦合干涉现象可以通过增大$\mathrm{d}n$来消除。这里要特别提醒不要为了追求完美的谱形而盲目加大$\mathrm{d}n$因为过大的耦合强度会让谐振峰分裂成多个峰更难处理。5.3 程序运行太慢怎么办如果波长扫描点设到5000个以上加上多层参数循环MATLAB跑起来确实会卡。我的经验是三步优化一是预分配数组不要在循环里动态扩展二是把循环内不变的量提到循环外比如$\kappa$的系数、$\Delta n_{eff}$的值三是考虑用向量化计算替代for循环把波长范围直接代入矩阵运算。% 向量化写法示例 lambda_vec lambda_range; kappa_vec pi * dn ./ lambda_vec * 0.8; delta_vec pi * (1/Lambda delta_n / lambda_vec); gamma_vec sqrt(kappa_vec.^2 delta_vec.^2); T_vec cos(gamma_vec * L).^2 (delta_vec.^2 ./ gamma_vec.^2) .* sin(gamma_vec * L).^2;这段向量化代码跑2001个波长点耗时从原来的约2秒降到不足0.1秒。虽然原始数据量不大、差距不明显但当你扩展到10000个点、还要嵌套周期和$\mathrm{d}n$扫描的时候这个优化能帮你节约大量时间。5.4 多峰分裂现象的排查当你发现谐振峰旁边出现两个甚至三个紧挨着的极小值时不用慌。这大概率是选定的包层模式不止一个。LPG和FBG的差异就在这里FBG只涉及纤芯基模的反向耦合模式单一LPG的包层模式有很多阶LP$_0$、LP$_1$、LP$_2$…… 每一阶都有自己的有效折射率和谐振波长。如果代码里没有明确限定只算第一个包层模式就可能同时算出好几个谐振峰。解决方法是给模式阶数加参数比如m 1然后只取对应阶数的$n_{eff}^{cl}$参与计算。如果你确实想仿真多个模式的叠加效果那就把每阶模式的透射谱分别算出来再按功率加权叠加而不是让代码在同一个循环里懵懵懂懂地混算。最后再分享一个小经验这个仿真代码做完之后不要急着删。把参数扫描的过程整理成一份简单的说明文档记录哪组参数对应哪条曲线后面写论文、补实验的时候会非常感谢自己当时做了这个记录。我自己就是这样几乎所有LPG相关的实验数据最后都能跟仿真曲线对得上。本文还有配套的精品资源点击获取