基于Matlab的光纤光栅仿真：从FBG/LPFG原理到传感器设计实践

简介：本资源是一套面向光学工程、光纤传感及光通信方向本科生与研究生的MATLAB仿真代码集，聚焦布拉格光纤光栅（FBG）与长周期光纤光栅（LPFG）的光谱特性建模与参数分析。资源解决初学者在理解耦合模理论、反射/传输谱计算及温度/应变响应机制时缺乏可运行实例的问题，覆盖均匀FBG、高斯型啁啾FBG、LPFG模式耦合等典型模型。压缩包共9个文件，全部为.m脚本，包含核心仿真函数（如fai_delta_Neff、gauss_function）、主调用程序（Uniform_FBG、CFBG_GaussianVisibility）及温度扰动分析模块（TemperatureChange系列），总大小仅7KB，轻量易读、结构清晰，便于逐行调试与参数修改。已有293人学习下载，读者可直接运行获取FBG尖锐反射峰、LPFG宽带损耗带等关键光谱图，快速掌握光栅周期、折射率调制深度、温度变化对中心波长漂移的影响规律，并复现文献中典型仿真结果。

1. 项目概述：从“光栅”到“传感”的桥梁

如果你正在研究光纤传感、光通信或者光学测量，那么“FBG”和“LPFG”这两个词一定不会陌生。它们全称分别是光纤布拉格光栅和长周期光纤光栅，是刻写在光纤纤芯里的周期性微结构，堪称光纤领域的“魔法刻痕”。我最初接触它们时，也被各种反射谱、透射谱、耦合模理论搞得头大，但当你真正上手用Matlab把它们“仿真”出来，亲眼看到理论公式如何变成屏幕上那条起伏的光谱曲线时，一切就豁然开朗了。这个项目，就是带你用Matlab这把“数字刻刀”，亲手“雕刻”出FBG和LPFG，并深入理解它们背后的物理机制、设计参数以及核心应用场景。无论你是光学工程的学生，还是初入光纤传感行业的工程师，这篇从理论推导到代码实现的完整指南，都能让你不仅知其然，更知其所以然，快速搭建起属于自己的光栅仿真与分析平台。

简单来说，FBG像一个波长特定的“镜子”，只反射某个很窄波段的光，对环境温度、应变极其敏感，是高性能点式传感的明星。而LPFG则像一个“波长过滤器”，将特定波长的光从纤芯耦合到包层并损耗掉，其透射谱对环境折射率变化非常灵敏，常用于生化传感。用Matlab仿真它们，核心目的有三个：一是在投入昂贵的实际刻写（如紫外激光曝光）前，进行充分的理论验证和参数优化，节省大量试错成本；二是深入理解光栅参数（如周期、长度、调制深度）如何影响最终的光谱特性，这是设计传感器的基本功；三是为后续的信号处理、解调算法开发提供一个可靠的、可任意调整的“数字孪生”测试环境。接下来，我们就从最根本的原理开始，一步步拆解如何用Matlab实现这两种光栅的仿真。

2. 核心原理与模型构建：光栅是如何“工作”的？

在动手写代码之前，我们必须先弄清楚FBG和LPFG到底是怎么改变光信号的。这离不开耦合模理论这个强大的工具。你可以把它理解为描述光如何在两个“通道”之间交换能量的数学语言。对于光栅，这两个“通道”就是光纤中传播的两种模式。

2.1 光纤布拉格光栅：纤芯模式的“自我对话”

FBG的原理相对直接。它通过紫外激光等手段，在光纤纤芯内引入一个周期性的折射率调制。这个周期Λ通常与光波长量级相当（例如~0.5微米）。当宽带光在光纤中传播时，满足布拉格条件的光会被强烈反射。

布拉格条件是理解FBG的钥匙：λ_B = 2 * n_eff * Λ。这里，λ_B就是布拉格波长（反射峰的中心波长），n_eff是光纤纤芯基模的有效折射率。这个公式意味着，FBG就像一个针对特定波长λ_B的谐振腔，只有波长满足这个条件的光，才会因为前后周期的反射光相干叠加而增强，最终被反射回来。

在耦合模理论框架下，FBG中发生的是反向耦合：纤芯中向前传播的基模（模式1）被光栅耦合成了向后传播的基模（模式2）。描述这一过程的是一组微分方程。最常用且高效的解法是传输矩阵法。它的核心思想是将整个光栅（长度L）离散成N个足够短的小段，每一小段可以视为一个均匀光栅。计算每一个小段的传输矩阵（一个2x2的矩阵，联系该段输入和输出的前向、后向光场），然后将所有小段的矩阵按顺序连乘，就能得到整个光栅的总体传输矩阵，进而求出反射谱和透射谱。

注意：为什么不用直接积分微分方程？因为传输矩阵法数值上更稳定、计算效率更高，尤其适合处理非均匀（如切趾、啁啾）光栅。它是工业界和学术界仿真FBG的事实标准。

2.2 长周期光纤光栅：纤芯到包层的“能量转移”

LPFG的周期Λ比FBG大得多，通常在几百微米量级。它导致的光学现象也完全不同：发生的是同向耦合。纤芯中向前传播的基模（LP01模）被耦合到同向传播的包层模（如LP11, LP12等）中。由于包层模很快被涂覆层损耗掉，所以在透射谱上会看到一个或多个损耗峰。

其相位匹配条件为：λ_res = (n_co - n_cl) * Λ。这里，λ_res是谐振波长（损耗峰位置），n_co和n_cl分别是纤芯基模和某个包层模的有效折射率。由于包层模的有效折射率对包层外介质（环境）的折射率非常敏感，所以λ_res会随着环境变化而漂移，这就是LPFG用于生化传感的物理基础。

仿真LPFG的模型比FBG稍复杂，因为涉及多个模式（一个纤芯模和多个包层模）之间的耦合。通常采用多模耦合模理论，并同样使用传输矩阵法求解。关键在于计算不同波长下，各模式之间的耦合系数。

2.3 仿真模型的关键参数与Matlab实现思路

无论是FBG还是LPFG，我们的Matlab模型都需要以下核心参数：

• 光栅结构参数：周期Λ、长度L、折射率调制深度Δn。这是设计的“输入”。
• 光纤参数：纤芯/包层折射率、半径，用于计算模式有效折射率n_eff。对于LPFG，还需要计算包层模的有效折射率。
• 仿真波段：起始波长λ_start、终止波长λ_end、波长点数。
• 输出结果：反射谱R(λ)（FBG）或透射谱T(λ)（LPFG）。

在Matlab中，我们将构建一个工作流：1) 设置参数；2) 计算波长相关的模式有效折射率（可能需要调用模式求解器或近似公式）；3) 对于每个波长，构建并连乘传输矩阵；4) 计算并绘制光谱。

3. Matlab仿真实现：从零搭建FBG仿真器

理论铺垫完毕，我们进入实战环节。我将以FBG为例，展示最核心的传输矩阵法代码实现。LPFG的代码结构类似，但模式计算部分更复杂。

3.1 环境准备与参数定义

首先，我们在Matlab脚本中定义所有基本参数。清晰的参数定义是后续调试的基础。

% FBG仿真参数设置 clear; close all; clc; % 1. 光纤参数 n_core = 1.4682; % 纤芯折射率 (SMF-28 @1550nm) n_clad = 1.4628; % 包层折射率 % 有效折射率近似，更精确的值应从模式求解获得 n_eff = 1.465; % 纤芯基模有效折射率近似值 % 2. 光栅参数 Lambda = 530e-9; % 光栅周期 (单位: 米)。根据布拉格公式，目标波长~1550nm L = 0.01; % 光栅长度 (单位: 米)，即1厘米 delta_n = 1e-4; % 折射率调制深度 (峰-峰值) % 3. 仿真光谱范围 lambda_start = 1540e-9; % 起始波长 (米) lambda_end = 1560e-9; % 终止波长 (米) num_points = 2000; % 波长点数 lambda = linspace(lambda_start, lambda_end, num_points)'; % 4. 传输矩阵法离散参数 N_segments = 200; % 将光栅离散成的段数 dz = L / N_segments; % 每段长度

实操心得：delta_n（折射率调制深度）是一个非常关键但常被初学者低估的参数。它直接影响光栅的反射率强度和带宽。对于普通传感FBG，1e-4量级是典型的；对于高反射率滤波器，可能需要1e-3或更高。一开始可以用1e-4作为基准进行调试。

3.2 传输矩阵法的核心代码实现

这是整个仿真的心脏。我们为每一个波长点计算整个光栅的响应。

% 初始化结果数组 R = zeros(size(lambda)); % 反射谱 T = zeros(size(lambda)); % 透射谱 % 主循环：遍历每个波长点 for i = 1:length(lambda) lam = lambda(i); % 计算该波长下的失谐量 (detuning) % 布拉格条件: lam_Bragg = 2 * n_eff * Lambda % 失谐量 delta = 2*pi*n_eff/lam - pi/Lambda % 更通用的形式： delta = (2*pi*n_eff/lam) - (pi/Lambda); % 计算耦合系数 kappa % 对于均匀正弦调制， kappa = pi * delta_n / lam kappa = pi * delta_n / lam; % 计算传播常数 gamma = sqrt(kappa^2 - delta^2) % 注意：当 kappa^2 < delta^2 时，gamma为虚数，对应带外区域 gamma_sq = kappa^2 - delta.^2; if gamma_sq >= 0 gamma = sqrt(gamma_sq); else gamma = 1i * sqrt(-gamma_sq); % 虚数单位 end % 初始化整体传输矩阵为单位矩阵 F_total = eye(2); % F = [F11, F12; F21, F22] % 遍历每一小段光栅 for seg = 1:N_segments % 计算单段均匀光栅的传输矩阵 (解析解) % 假设这一段内参数恒定 s = gamma * dz; if gamma_sq >= 0 % 带内 (通带) F11 = cosh(s) - 1i*(delta/gamma)*sinh(s); F12 = -1i*(kappa/gamma)*sinh(s); F21 = 1i*(kappa/gamma)*sinh(s); F22 = cosh(s) + 1i*(delta/gamma)*sinh(s); else % 带外 (阻带) % 此时 gamma 为虚数，令 gamma = 1i*|gamma|，代入公式会变为正弦余弦形式 gamma_abs = abs(gamma); F11 = cos(gamma_abs*dz) - (delta/gamma_abs)*sin(gamma_abs*dz); F12 = -(kappa/gamma_abs)*sin(gamma_abs*dz); F21 = (kappa/gamma_abs)*sin(gamma_abs*dz); F22 = cos(gamma_abs*dz) + (delta/gamma_abs)*sin(gamma_abs*dz); end F_seg = [F11, F12; F21, F22]; % 连乘得到整体矩阵 F_total = F_seg * F_total; end % 从整体传输矩阵计算反射系数和透射系数 % 边界条件：光从左侧入射，右侧无入射光 r = -F_total(2,1) / F_total(2,2); % 反射系数 (复数) t = F_total(1,1) + F_total(1,2)*r; % 透射系数 (复数) % 计算反射率和透射率 (光强) R(i) = abs(r).^2; T(i) = abs(t).^2; end

3.3 结果可视化与基础分析

有了反射谱R和透射谱T，我们可以直观地看到FBG的特性。

% 绘制反射谱和透射谱 figure('Position', [100, 100, 1200, 500]); subplot(1,2,1); plot(lambda*1e9, R*100, 'b-', 'LineWidth', 1.5); xlabel('波长 (nm)'); ylabel('反射率 (%)'); title('FBG反射光谱'); grid on; xlim([lambda_start, lambda_end]*1e9); % 标记布拉格波长 lambda_Bragg = 2 * n_eff * Lambda; hold on; plot([1,1]*lambda_Bragg*1e9, [0, max(R*100)], 'r--', 'LineWidth', 1); legend('反射谱', sprintf('布拉格波长 %.2f nm', lambda_Bragg*1e9)); subplot(1,2,2); plot(lambda*1e9, T*100, 'r-', 'LineWidth', 1.5); xlabel('波长 (nm)'); ylabel('透射率 (%)'); title('FBG透射光谱'); grid on; xlim([lambda_start, lambda_end]*1e9); hold on; plot([1,1]*lambda_Bragg*1e9, [0, 100], 'b--', 'LineWidth', 1); % 计算并显示关键光谱参数 [maxR, idx] = max(R); lambda_max = lambda(idx) * 1e9; fwhm = find(lambda >= lambda_max*1e-9/2 & R >= maxR/2); % 简化的半高全宽查找 if length(fwhm) >= 2 bandwidth = (lambda(fwhm(end)) - lambda(fwhm(1))) * 1e9; else bandwidth = NaN; end fprintf('===== FBG仿真结果 =====\n'); fprintf('理论布拉格波长: %.3f nm\n', lambda_Bragg*1e9); fprintf('仿真反射峰中心波长: %.3f nm\n', lambda_max); fprintf('峰值反射率: %.2f%%\n', maxR*100); fprintf('反射谱-3dB带宽: %.3f nm\n', bandwidth);

运行这段代码，你将得到一张清晰的反射/透射光谱图，并能在命令行看到计算出的关键参数。这构成了FBG仿真的最基础、最核心的框架。

4. 进阶仿真：非均匀FBG与LPFG初步

掌握了均匀FBG的仿真，我们就可以探索更复杂、更实用的结构。

4.1 实现切趾与啁啾FBG

实际应用中，为了抑制反射谱的旁瓣，我们会使用切趾，即让折射率调制深度Δn沿光栅长度方向按一定函数（如高斯、升余弦）变化。为了得到宽带的反射谱或用于色散补偿，会使用啁啾，即让光栅周期Λ沿长度方向线性或非线性变化。

在代码中实现这些特性非常直观：只需在分段循环内，根据当前段的位置z，动态计算该段的delta_n或Lambda即可。

% 示例：高斯切趾 apodization_profile = exp(-(( (1:N_segments)-N_segments/2) / (N_segments/4)).^2); % 在分段循环内： delta_n_local = delta_n * apodization_profile(seg); % 示例：线性啁啾 chirp_rate = 1e-9; % 啁啾系数，单位：米/米 (周期随长度的变化率) % 在分段循环内： Lambda_local = Lambda + chirp_rate * (seg - 0.5) * dz; % 当前段的周期 % 同时需要重新计算当前段的失谐量 delta_local

4.2 LPFG仿真框架搭建

LPFG的仿真框架与FBG类似，但核心区别在于模式计算。你需要计算一系列包层模的有效折射率n_clad,m。这通常需要求解光纤的标量或矢量波动方程，对于标准单模光纤，有经验公式或查表可用。为了简化演示，我们假设已知一个包层模的有效折射率n_clad。

% LPFG 参数设置 Lambda_lpfg = 400e-6; % LPFG周期，单位：米，典型值几百微米 L_lpfg = 0.03; % 光栅长度，通常比FBG长 delta_n_lpfg = 5e-4; % 调制深度 % 假设的包层模有效折射率 (需要根据实际模式计算) n_clad_mode = 1.45; % 示例值 % 初始化LPFG透射谱 T_lpfg = ones(size(lambda)); for i = 1:length(lambda) lam = lambda(i); % 计算纤芯模和包层模的传播常数 beta_co = 2*pi*n_eff / lam; beta_cl = 2*pi*n_clad_mode / lam; % LPFG的相位失谐量 (同向耦合) delta_lpfg = (beta_co - beta_cl) / 2 - pi / Lambda_lpfg; % 耦合系数 (近似) kappa_lpfg = pi * delta_n_lpfg / lam; % 注意：实际耦合系数计算更复杂 % 使用传输矩阵法 (公式与FBG类似，但耦合机制不同，此处简化) % 对于均匀LPFG，透射率有解析解: T = 1 - [kappa^2 / (kappa^2 + delta^2)] * sin^2( sqrt(kappa^2+delta^2) * L ) gamma_lpfg_sq = kappa_lpfg^2 + delta_lpfg^2; if gamma_lpfg_sq > 0 gamma_lpfg = sqrt(gamma_lpfg_sq); T_lpfg(i) = 1 - (kappa_lpfg^2 / gamma_lpfg_sq) * (sin(gamma_lpfg * L_lpfg))^2; end end % 绘制LPFG透射谱 figure; plot(lambda*1e9, T_lpfg*100, 'LineWidth', 1.5); xlabel('波长 (nm)'); ylabel('透射率 (%)'); title('LPFG透射光谱 (单模耦合简化模型)'); grid on; xlim([lambda_start, lambda_end]*1e9);

这个LPFG模型是高度简化的。一个完整的LPFG仿真器需要：1) 准确的包层模有效折射率计算（可能涉及贝塞尔函数求根）；2) 更精确的耦合系数模型（与光栅类型和模式重叠积分有关）。这通常需要更专业的工具箱或更复杂的代码。

5. 仿真结果分析与传感器设计启示

通过运行上述代码，我们得到了光谱。如何解读这些曲线，并用于指导传感器设计呢？

5.1 FBG光谱参数解读与传感原理

对于FBG，我们最关心反射谱的中心波长λ_B、峰值反射率R_max和带宽FWHM。

• 传感原理：当FBG受到轴向应变ε或温度变化ΔT时，其周期Λ和有效折射率n_eff会发生变化，导致λ_B漂移。灵敏度公式为：Δλ_B / λ_B = (1 - p_e) * ε + (α + ξ) * ΔT其中p_e是弹光系数，α是热膨胀系数，ξ是热光系数。对于1550nm附近的石英光纤FBG，典型的应变灵敏度约为~1.2 pm/με，温度灵敏度约为~10 pm/°C。
• 设计启示：
  • 反射率：由κL乘积决定（κ是耦合系数，L是长度）。κL越大，反射率越高。要达到高反射率（>90%），需要较大的Δn或较长的L。
  • 带宽：近似公式为Δλ ≈ (λ_B^2)/(π * n_eff * L) * sqrt((κL)^2 + π^2)。带宽越窄，波长检测分辨率要求越高，但传感器串扰越小。
  • 旁瓣：均匀FBG旁瓣很高，可能造成传感解调时的串扰。必须使用切趾来抑制旁瓣。高斯切趾和升余弦切趾是常用选择。

5.2 LPFG光谱参数解读与传感原理

对于LPFG，我们关注透射谱中损耗峰的中心波长λ_res、深度和**宽度`。

• 传感原理：LPFG的谐振波长λ_res主要取决于纤芯模与包层模的有效折射率差(n_co - n_cl)。包层模的场分布在包层中，其有效折射率对外界包层材料的折射率变化极其敏感。因此，LPFG是折射率传感器的天然选择。此外，温度、应变也会通过改变光纤本身性质影响λ_res，但通常折射率灵敏度远高于温漂。
• 设计启示：
  • 周期选择：根据目标传感波长和模式，由相位匹配条件Λ = λ_res / (n_co - n_cl)初步确定周期。
  • 耦合强度：损耗峰深度由耦合系数κ和光栅长度L决定。要达到深度耦合（接近100%损耗），同样需要优化κL。
  • 模式选择：不同阶次的包层模（LP11, LP12, ...）具有不同的场分布和对外界折射率的灵敏度。高阶包层模通常具有更高的灵敏度，但其谐振峰可能更宽。

5.3 利用仿真进行传感器优化

Matlab仿真的真正威力在于参数扫描和优化。你可以轻松地写一个循环，研究某个参数（如L、Δn）对光谱的影响。

% 示例：研究光栅长度L对FBG反射谱的影响 lengths = [0.005, 0.01, 0.02, 0.04]; % 单位：米 figure; hold on; for L_val = lengths % 调用之前写好的FBG仿真函数，传入新的L_val % 假设你已将FBG仿真封装为函数 [R, lambda] = simulate_FBG(n_eff, Lambda, L_val, delta_n, ...) [R_temp, ~] = simulate_FBG(n_eff, Lambda, L_val, delta_n, lambda); plot(lambda*1e9, R_temp*100, 'DisplayName', sprintf('L=%.1f cm', L_val*100)); end xlabel('波长 (nm)'); ylabel('反射率 (%)'); title('不同光栅长度下的FBG反射谱'); legend('show'); grid on;

通过这样的分析，你可以在制作实物前，就确定满足你传感需求（如特定反射率、带宽）的最佳光栅参数组合，极大提升研发效率。

6. 常见问题、调试技巧与高级扩展

在实际仿真和后续应用中，你肯定会遇到各种问题。这里分享一些我踩过的坑和解决思路。

6.1 仿真结果异常排查表

6.2 性能优化与代码进阶

• 向量化与预计算：尽可能将循环内的计算向量化。例如，可以一次性计算所有波长点的delta和kappa向量，然后在分段循环外进行部分矩阵运算。
• 函数封装：将FBG和LPFG的仿真核心代码封装成函数，例如[R, T, lambda] = simulateUniformFBG(params)。这样主脚本会非常清晰，也便于参数化扫描。
• 引入模式求解器：要获得更精确的结果，尤其是对于LPFG或特殊光纤，需要集成模式求解。可以自己编写基于标量近似或矢量解的模式求解代码，或者利用一些开源工具箱。
• GUI开发：使用Matlab的App Designer，可以快速搭建一个带有滑块、输入框和实时绘图的光栅仿真GUI，这对于教学和快速原型演示非常有用。

6.3 从仿真到实际应用的联系

仿真最终要服务于实践。当你设计一个FBG温度传感器时，仿真可以告诉你：为了在1550nm附近获得80%的反射率和0.2nm的带宽，你需要多长的光栅和多大的折射率调制。然后，你可以将这些参数反馈给光栅刻写工艺（如相位掩模法）的工程师。

对于LPFG生化传感器，仿真可以帮助你预测：当包层外介质折射率从1.33（水）变化到1.45（某种蛋白质溶液）时，损耗峰会漂移多少纳米。这直接决定了你后续解调系统需要多大的波长检测范围和多高的分辨率。

最后，这个仿真平台还可以作为解调算法的测试床。你可以轻易地在仿真出的理想光谱上添加高斯白噪声来模拟探测器噪声，或者模拟多个FBG的谱叠加来测试波分复用/时分复用解调算法的性能。

本文还有配套的精品资源，点击获取